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ABSTRACT 

We use a disk galaxy evolution model to investigate the impact of mass outflows (a.k.a. 
feedback) on disk galaxy scaling relations, mass fractions and spin parameters. Our 
model follows the accretion, cooling, star formation and ejection of baryonic mass 
inside growing dark matter haloes, with cosmologically motivated angular momentum 
distributions, and dark matter halo structure. Models without feedback produce disks 
that are too small, too gas poor, and which rotate too fast. Feedback reduces the galaxy 
formation efficiency, cgf, (defined as the fraction of the universally available baryons 
that end up as stars and cold gas in a given galaxy) , resulting in larger disks with higher 
gas fractions and lower rotation velocities. Models with feedback can reproduce the 
zero points of the scaling relations between rotation velocity, stellar mass and disk size, 
but only in the absence of adiabatic contraction. Our feedback mechanism is maximally 
efficient in expelling mass, but our successful models require 25% of the SN energy, or 
100% of the SN momentum, to drive an outflow. It remains to be seen whether such 
high efficiencies are realistic or not. Our energy and momentum driven wind models 
result in different slopes of various scaling relations. Energy driven winds result in 
steeper slopes to the galaxy mass - halo mass, and stellar mass - halo mass relations, a 
shallower slope to the galaxy size - stellar mass relation at z = 0, and a steeper slope to 
the cold gas metallicity - stellar mass relation at z ~ 2. Observations favor the energy 
driven wind at stellar masses below M stal < 10 10 5 Mq, but the momentum driven wind 
model at high masses. The ratio between the specific angular momentum of the baryons 
to that of the halo, (jgai/w ga i), is not unity in our models with inflow and outflow. 
Yet this is the standard assumption in models of disk formation. Above a halo mass of 
M v ; r ~ 10 12 A/q, cooling becomes increasingly inefficient, which results in (j ga i/m ga i) 
decreasing with increasing halo mass. Below a halo mass of M v - lr ~ 10 12 M Q , feedback 
becomes increasingly efficient. Feedback preferentially ejects low angular momentum 
material because star formation is more efficient at smaller galactic radii, and at higher 
redshifts. This results in (jgai/'rigai) increasing with decreasing halo mass. This effect 
helps to resolve the discrepancy between the high spin parameters observed for dwarf 
galaxies with the low spin parameters predicted from ACDM . 
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1 INTRODUCTION 

Understanding the origin and nature of galaxy scaling re- 
lations is a fundamental quest of any successful theory of 
galaxy formation. The success of a particular theory will be 
judged by its ability to reproduce the slope, scatter, and 
zero-point of any robust galaxy scaling relation. Of particu- 
lar interest are the scaling relations between size (R), lumi- 
nosity (or stellar mass, M) (L), and velocity (V), as these 
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parameters are related to each other via the virial theorem. 
For early type galaxies, these three parameters are indeed 
coupled, resulting in a two-dimensional plane known as the 
Fundamental Plane (FP: Dressier et al. 1987; Djorgovski & 
Davis 1987). However, for late-type galaxies, the relation be- 
tween rotation velocity and luminosity, known as the Tully- 
Fisher relation (TF: Tully & Fisher 1977) is independent of 
galaxy size or surface brightness (e.g. Courteau & Rix 1999; 
Courteau et al. 2007; Pizagno et al. 2007). Understanding 
the origin of this surface brightness independence is likely 
the key to understanding the small scatter of the TF, and 
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may even explain the origin of the zero point of the TF 
relation. 

The slopes of the VLR relations for disk galaxies can 
be broadly understood with galaxy formation models that 
include virial equilibrium after dissipation-less collapse of 
quasi-spherical cold dark matter (CDM) halos and angular 
momentum conservation (e.g. Mo, Mao, & White 1998; van 
den Bosch 1998, 2000; Navarro & Steinmetz 2000; Firmani 
& Avila-Reese 2000; Dutton et al. 2007). 

However, reproducing these scaling relations in detail 
has been a problem for galaxy formation models. To date, 
no (self-consistent) CDM-based model of galaxy formation 
can simultaneously match the the zero points of the TF rela- 
tion, galaxy sizes, and the luminosity (or stellar mass) func- 
tion (e.g. Benson et al. 2003; Dutton et al. 2007). Models 
which match the zero point of the TF relation, do so by 
making the assumption that, V TO t, the observed rotation ve- 
locity is equal to, V v i r , the circular velocity at the virial ra- 
dius (e.g. Somerville & Primack 1999), or, V max ,h, the maxi- 
mum circular velocity of the halo prior to galaxy formation, 
(e.g. Croton et al. 2006). For typical galaxy mass dark mat- 
ter haloes Knax.h/K/ir — 1-1 (e.g. Bullock et al. 2001a), so 
these two assumptions are almost equivalent. Other more 
observational approaches also support the conclusion that 
Kot ^ Knax.h (Eke et al. 2006; Blanton et al. 2008). The 
problem for galaxy formation theory is that both cosmolog- 
ical simulations and analytic models of disk galaxy forma- 
tion that take into account the self-gravity of the baryons 
and the effect of halo contraction (Blumenthal et al. 1986) 
find that V TO t ^ 1.8V v ir (Navarro & Steinmetz 2000; Abadi 
et al. 2003; Dutton et al. 2007; Governato et al. 2007). 

Due to the almost one-to-one correlation between 
Vmax.h and M v i r , and the correlation between Vrot and bary- 
onic Mgai mass (the baryonic TF relation e.g. McGaugh 
et al. 2000) an equivalent constraint to the ratio between Vrot 
and Vmax.h is the galaxy mass fraction: m ga i = Af ga i/M v i r . 
Observations of halo masses using weak lensing (Hoek- 
stra et al. 2005; Mandelbaum et al. 2006) find that for 
late type galaxies the maximum m ga i — 0.33/bar, where 
/bar = fib/^m.o — 0.17 is the cosmic baryon mass fraction 
at redshift zero. Similarly low galaxy formation efficiencies 
are obtained by methods that match the stellar mass func- 
tion with the dark halo mass function (e.g. Yang et al. 2007; 
Conroy & Wechsler 2008). 

A partial explanation for the surface brightness inde- 
pendence of the TF relation, or equivalently the weak corre- 
lation between residuals of the VL and RL relations, is that 
observationally (e.g. McGaugh & de Blok 1997) gas frac- 
tions correlate with surface brightness, with higher gas frac- 
tions in lower surface brightness galaxies. Since lower surface 
density disks are expected to rotate more slowly at a given 
baryonic mass, the larger gas fractions shift these galaxies to 
lower stellar masses and hence lower luminosities (Firmani 
& Avila-Reese 2000; van den Bosch 2000). However, this so- 
lution is not very effective for the high surface brightness, 
low gas fraction galaxies. Dutton et al. (2007) showed that 
a reasonable amount of scatter in the stellar mass-to-light 
ratios helps to reduce the correlation between the residuals 
of the VL and RL relations further. Gnedin et al. (2007) 
proposed a correlation between disk mass fraction and disk 
surface density to explain the lack of a correlation between 
the VM and RM relations. 



Thus understanding the physical mechanisms that de- 
termine galaxy mass fractions is fundamental to our under- 
standing of the origin of the VLR relations of disk galaxies. 
In the standard picture of disk galaxy formation, gas that 
enters the halo gets shock heated to the virial temperature, 
and then cools radiatively (Fall & Efstathiou 1980). Thus in 
order to produce low galaxy mass fractions, either a signif- 
icant fraction of the gas has to be prevented from cooling, 
or a significant fraction must subsequently be ejected from 
the disk and halo. The latest hydrodynamical simulations 
indicate that rather than accreting gas via a cooling flow, 
below a critical halo mass of M v i r ~ 1O 12 M0 gas does not 
shock as it enters the halo, and instead accretes straight 
onto the galaxy in a cold flow (Birnboim & Dekel 2003; 
Keres et al. 2005). In this scenario, essentially all of the 
baryons accrete onto the galaxy, and the problem of stop- 
ping the baryons from cooling becomes one of stopping the 
cold flows from forming. Thus in the absence of some kind 
of pre-heating (e.g. Mo et al. 2005) which would shut down 
the cold flows, mass outflows are required in order to pro- 
duce the low galaxy mass fractions observed in the Universe 
today. 

In order to investigate how mass outflows (a.k.a. feed- 
back) determines galaxy mass fractions, and the impact this 
has on disk galaxy scaling relations, we use an updated ver- 
sion of the disk galaxy evolution models presented in van den 
Bosch (2001; 2002a). In these models the input parameters 
are the concentration and spin of the dark matter halo. The 
galaxy spin parameter, galaxy mass fraction and gas to stel- 
lar mass ratio are collectively determined by the efficiencies 
of cooling, star formation, and feedback. 

The main differences with respect to the van den Bosch 
(2002a) models are: 

(i) We use cosmologically motivated specific angular mo- 
mentum distributions (AMDs) of the halo gas rather than 
shells in solid body rotation; 

(ii) We consider a star formation recipe based on dense 
molecular gas, rather than on total gas with a Toomre star 
formation threshold; 

(iii) We include scatter in halo concentration which we 
relate to the mass accretion history (MAH); 

(iv) We explore two different feedback models: one based 
on kinetic energy conservation, the other based on momen- 
tum conservation. 

An important aspect of this model is that we do not 
assume that the baryonic disk has an exponential density 
profile. In this model the surface density profile of the bary- 
onic disk is determined by detailed conservation of angular 
momentum, starting from the AMDs of gas haloes as found 
in cosmological simulations. The surface density profile of 
the stars is then determined by the relative efficiencies of 
star-formation, outflows, and inflows as a function of ra- 
dius. This allows us to self-consistently follow the evolution 
of the radial distributions of gas and stars. In future papers 
we use this galaxy formation model to investigate the origin 
of disk galaxy surface density profiles (Dutton 2008) and the 
evolution of disk galaxy scaling relations. 

This paper is organized as follows: In Sj2]we describe the 
disk galaxy formation models; in fj3] we discuss the effect 
of feedback on the VMR relations; in i|4] we discuss the 
impact of feedback on galaxy mass fractions and galaxy spin 
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parameters; in Sj5] we discuss the effect of feedback on the 
mass-metallicity relation at z ~ 2; in f|6] we discuss how 
galaxies lose their gas; and in fJ7]we give a summary. 



2 DISK GALAXY FORMATION MODELS 

The main assumptions that characterize the framework of 
our models are the following: 

(i) Dark matter haloes around disk galaxies grow by the 
smooth accretion of mass; 

(ii) The baryons acquire the same distribution of specific 
angular momentum as the dark matter; 

(iii) Gas that enters the halo is shock heated to the virial 
temperature; 

(iv) Gas cools conserving its specific angular momentum; 

(v) Star formation occurs according to a Schmidt type 
law on the dense molecular gas; 

(vi) Supernova feedback re-heats some of the cooled gas, 
ejecting it from the halo; 

(vii) The halo contracts and expands adiabatically to in- 
flows and outflows, respectively; 

Each model galaxy consists of five mass components: 
dark matter, hot halo gas, disk mass in stars, disk mass 
in cold gas, and ejected gas. The dark matter and the hot 
gas are assumed to be distributed in spherical shells, the 
disk mass is assumed to be in infinitesimally thin annuli. 
Throughout this paper we refer to R as radius, t as time 
(where t — is defined as the Big Bang) and z as redshift. 

For each galaxy we set up a radial grid with 200 bins 
quadratically sampled from between 0.001 to 1 times the 
redshift zero virial radius, and we follow the formation and 
evolution of the galaxy using 400 redshift steps quadrati- 
cally sampled from z — 10 to z = 0. For each time step 
we compute the changes in the various mass components in 
each radial bin. The prescriptions we use are described in 
detail below. 

2.1 Limitations of the Model 

Before we describe the details of the model, we first discuss 
the relevance of our model to understanding galaxy forma- 
tion in a hierarchical universe. 

The assumption of smooth mass accretion might seem 
inconsistent with the hierarchical merger picture of struc- 
ture formation in a ACDM universe. However, major merg- 
ers of stellar rich galaxies are known to destroy disks (e.g. 
Barnes 1992; Cox et al. 2006), so disk galaxies are unlikely 
to form in haloes with recent major mergers. Minor mergers 
are more common than major mergers, and are likely to play 
and important role in the formation of galaxy bulges, either 
directly, or by triggering secular processes. Thus by not in- 
cluding mergers our model underestimates bulge fractions. 
However, one of the goals of these models is to determine 
how much of the structural properties of disk galaxies can 
be accounted for with the "zeroth order" scenario of disk 
formation (smooth accretion and quiescent star formation). 

Our assumption about the way gas is accreted into 
galaxies (by a cooling flow of shock heated gas) is likely in- 
correct. Simulations suggest that disk galaxies accrete most 
of their mass though cold flows, and that in the absence of 



extra heating or outflows the baryon fraction of galaxies is 
close to the universal value (Keres et al. 2005). In our model 
the gas shock heats, but since for haloes with masses below 
10 12 M© the cooling time is short compared to the Hubble 
time, essentially all the gas that enters the halo accretes 
onto the disk in a free fall time. Thus although the phys- 
ical mechanism by which galaxies accrete their gas in our 
model and cosmological simulations are different, we expect 
that the specific angular momentum distribution of cold flow 
baryons to be the same as that of the dark matter, though 
this needs to be verified using cosmological simulations. 

Our assumption about the efficiency at which gas is 
converted into stars is necessarily over-simplified, but it is 
an improvement over the majority of star formation recipes 
(which may be physically or empirically motivated) used 
in semi-analytic models and hydro-dynamical simulations, 
which are based on the density (or mass) of the total cold gas 
(atomic plus molecular). We assume an empirical relation 
between the local star formation rate and the local density 
of molecular gas. We calculate the molecular gas fraction (as 
a function of radius in the disk) using the empirical relation 
between mid plane pressure and molecular gas fraction from 
Blitz & Rosolowski (2006). A more realistic treatment of 
star formation would model the formation and destruction of 
molecular gas in a physically motivated way (e.g. Pelupessy, 
Papadopoulos, & van der Werf 2006; Robertson & Kravtsov 
2008) 

In our feedback model we only consider winds that are 
able to escape the halo, and we assume that mass in such 
winds is lost forever. In reality mass that escapes the halo 
may fall back at later times. Furthermore, there may be 
winds that have enough energy /momentum to escape the 
disk but not the halo. The gas in these winds could then re- 
cool back onto the disk producing a galactic fountain. Since 
very little is known about how feedback works, and our main 
interest is to determine how much mass can be ejected from 
the disk and halo, rather than introducing additional free 
parameters to our wind model, we assume the maximal mass 
loss model. 

The assumption that the halo responds adiabatically 
to inflows and outflows may not be correct. When galaxies 
accrete their gas via a smooth cooling flow the gas radiates 
away its energy. However, when a galaxy acquires its gas 
via cold flows, clumps of cold gas can exchange energy with 
the halo via dynamical friction (e.g. El Zant, Shlosman, & 
Hoffman 2001; El-Zant et al. 2004; Jardel & Sellwood 2009), 
causing the halo to expand. In both scenarios the natural 
response of the halo to the deepening of the potential well 
due to the condensation of baryons is to contract, but in the 
latter scenario, the transfer of energy between baryons and 
dark matter will counter this effect. Processes internal to 
disks, such as bars, can also transfer energy to the halo via 
dynamical friction, causing it to expand (Weinberg & Katz 
2002; Holley-Bockelmann, Weinberg & Katz 2005; Sellwood 
2008). Mass outflows can also result in halo expansion, for 
example an adiabatic inflow followed by an instantaneous 
outflow can result in net halo expansion (Navarro, Eke, & 
Frenk 1996; Gnedin & Zhao 2002; Read & Gilmore 2005). 
Thus in order to asses how much impact the assumption of 
adiabatic inflow and outflow has on the structural properties 
of the resulting galaxies, we also consider models in which 
the halo does not respond to galaxy formation. 
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2.2 Dark Matter Haloes 

The backbone of our galaxy formation model is the growth 
of the dark matter halo, which we model by a smooth mass 
accretion history (MAH). Van den Bosch (2002b) and Wech- 
sler et al. (2002) have shown that the MAH is essentially a 
one parameter family. The MAH from Wechsler et al. (2002) 
is given by 



M vir (z) = M vir e 



(1) 



where M V [ T is the redshift zero mass and a is related to the 
epoch of formation via 

a = a c S. (2) 

Here a c — (1 + z c ) _1 is defined as the scale factor a when 
the logarithmic slope of the accretion rate falls below some 
specified value, S. Following Wechsler et al. (2002) we adopt 
S = 2. Before we discuss how to compute a c we discuss how 
the structural properties of the halo depend on its mass. 

In the standard spherical top-hat collapse model the 
virial radius, R v i I (z), of a halo of mass M v i T (z) at a redshift, 
z, is given by 



h *kpc 
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where h = ifo/100, and A v i r is the virial density, relative to 
the critical density for closure. We use the fitting formula of 
Bryan & Norman (1998) 

Avir = 18tt 2 + 82x - 39x 2 , (4) 
1. The evolution of the matter density is 



with x = Slm(z) 
given by 



n m ( z ) = n m , (i + zf 



H{z) 



H 



(5) 



and the evolution of the Hubble parameter is given by 



H 2 {z) = Hi [fiA + (l-n A -a^,o)(l + zr +On,0(l + 2) d ].(6) 

The relation between the virial velocity, V v \ r , and virial ra- 
dius, i?vir is given by 



GM vir 
-Rvir 



(7) 



We assume that the density profile of the halo at each 
redshift is given by an NFW (Navarro, Frenk, & White 1997) 
profile 



p(R) 



(8) 



Pent {R/R s )(l + R/R s ) 2 ' 

where R s is the radius where the slope of the density profile 
is -2, the so called scale radius, p cr it is the critical density of 
the universe, and 5 C is the characteristic overdensity of the 
halo. The concentration parameter of the halo is defined as 
c = Rvir/Rs and is related to <5 C via 

5c = A-tt: rr. (9) 



3 [ln(l + c) - c/(l + c)] ' 

Following Bullock et al. (2001a) we assume that the 
concentration parameter evolves as 

.1 + Zc 



c(z) = Ka{z)/a c = K 



l + z 



(10) 



with a minimum value of c(z) — K, corresponding to a con- 
stant concentration during the rapid accretion phase of dark 
halo growth (Zhao et al. 2003). Thus if we specify the red- 
shift zero concentration we can compute the collapse epoch 
via a c = K/c, and hence the MAH through Eqs. (0 & ©. 

To compute the mean concentration for a halo of a 
given mass at z = we use the Bullock et al. (2001a) 
model. This model requires us to specify the cosmology 
Qa, fim,o, ^b, °"8, n, h, as well as 2 free parameters, F and 
K. We assume the concentration is log-normally distributed 
with median c(M) from the Bullock et al. (2001a) model 
and scatter ci nc . Our adopted values for these parameters 
are given in £12.101 



2.3 Gas Cooling 

At each time step a shell with mass AM — M v i T (t) — M v i T (t — 
At) virializes. A fraction /bar = Slb/^m.o of this mass is in 
baryons, and is heated to the halo's virial temperature 



_ 1(0% 2 
vir ~2 k vir 



(11) 



where /im p is the mass per particle, and k is Boltzmann's 
(3) constant. The baryons dissipate energy radiatively, lose pres- 
sure support, and collapse until they reach centrifugal equi- 
librium. The time scale over which this occurs is given by 
t c = max[tf3F,t C ooi]- Here tg is the free fall time defined as 



tff = 



3tt 



(12) 



32Gp 

with p the average halo density, and 

, 3 kT v j r fl c /io\ 

tcooi = -Mm? — — — (13) 

£ PhotAN(,^hot J pc — 1 

is the cooling time. Here phot is the density of the hot gas, p e 
is the number of particles per electron, and An is the normal- 
ized cooling function for a gas with metallicity Zhot ■ For phot 
we use /bar times the density at the virial radius at the time 
the gas enters the halo. For An we use the collisional ioniza- 
tion equilibrium cooling functions of Sutherland & Dopita 
(1993), assuming a helium mass abundance of 0.25. 

For each time step, we compute the range of times be- 
tween which gas that collapses onto the disk in the cur- 
rent time interval, entered the halo. We label these times as 

^c,min(^) and £c,max(i). 



2.4 Angular Momentum Distribution 

The radius at which the cooled gas ends up depends on 
its specific angular momentum, j. Van den Bosch (2001; 
2002a) assumed the j-distribution to be that of a shell in 
solid body rotation. The angular momentum of this shell 
can be computed assuming that the spin parameter, A, is 
constant between time steps. The spin parameter is defined 
by 



A = 



7 • \E ■ l 1/2 

<J vir | -L^vir | 

gmX 2 



(14) 



where M v ir, Jvir, and i? v ir are the mass, total angular mo- 
mentum and energy of the halo, respectively. As shown by 
van den Bosch (2001) this results in density profiles that are 
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more concentrated than exponential in the center, and also 
that truncate at shorter radii than is observed in some disk 
galaxies. 

Building on the work of Bullock et al. (2001b), Sharma 
& Steinmetz (2005, hereafter SS05) used a series of non ra- 
diative N-body/SPH simulations in a ACDM cosmology to 
study the growth of angular momentum in galaxy systems. 
They introduced a function that is able to describe, with a 
single parameter, a, the specific angular momentum distri- 
bution of the gas and dark matter in their simulations, as 
well as that of exponential disks in NFW haloes: 



M(< j) 

Mvir 



7(Q, 



3d 



3d 



t/vir 1 



Mvir a 



(15) 



where 7 is the incomplete gamma function. SS05 found that 
the distribution of a is log-normally distributed with mean 
log a ~ —0.05 and standard deviation in logo ~ 0.11. 

In summary the distribution of specific angular momen- 
tum of the dark matter halo and hot gas can be described by 
two parameters: a normalization (A) and a shape (a). Both 
the normalization and shape parameters are log-normally 
distributed, with significant scatter. We assume that the 
spin and shape parameters are uncorrelated, although Bul- 
lock et al. (2001b) show that there may be a weak correlation 
between A and a. Furthermore we assume that, for a given 
halo, A and a are constant with redshift. These assumptions 
are made for simplicity, and need to be tested with cosmo- 
logical simulations. 

In order to compute the amount of mass, with a given 
specific angular momentum, j, that has collapsed onto the 
disk in each time step, M c , we take the difference between 
the distributions of specific angular momentum of the halos 
at times t c ,mm and t c ,max (t c ,min and t c ,max are defined in 
] above); 

-M vir (< j)(t c ,min)].(16) 



Mc(< j)(t) = /bar [Mvir (< j){tc,n 



2.5 Conservation of Angular Momentum and 
Halo Contraction 

In order to compute the radius, R, at which material with 
specific angular momentum, j, ends up we assume specific 
angular momentum is conserved, i.e. one should solve 



j =RV ci rc(R) 



(17) 



for R. Here Kirc is the total circular velocity (from stars, 
cold gas, hot gas, and dark matter). 

As the galaxy grows over time the circular velocity at 
a given radius increases. Thus to conserve specific angular 
momentum, material that settled at radius R would need to 
drift to smaller radii over time. Given that the gas and stars 
effect the circular velocity, strictly conserving the specific 
angular momentum of the baryons over time is difficult to 
implement numerically. To get around this problem we use 
specific angular momentum, j, rather than radius as our 
radial coordinate. 

Under the simplifying assumption that Kirc = [GM(< 
R)/R] 1/2 , where M(< R) is the total mass within a spherical 
radius, R, the radius that corresponds to a given j is given 
by 



R 



3 



GM{< j) 



(18) 



This has a number of desirable properties: 1) At each time 
step it is trivial to calculate how much cold gas is added to 
each bin in j. 2) Over time, as the potential well changes, the 
specific angular momentum of the baryons is automatically 
conserved; 3) The response of the halo to the cooling of the 
baryons is automatically taken into account. 4) The resulting 
radial grid is adaptive, as the mapping between j and R 
depends on the amount of mass enclosed. 

The price that we pay for these advantages is that the 
circular velocity due to the disk is not calculated correctly. 
Due to the disk geometry the true circular velocity deviates 
from that given by the spherically enclosed mass. For ex- 
ample the peak V c ii C of a thick exponential disk is ~ 10% 
higher than that obtained using the enclosed mass, and at 
small radii the proper V c irc increases linearly with radius, 
whereas using the enclosed mass V c i IC scales as 7? 1//2 . 

However, given that computing the proper circular ve- 
locity of the disk is very time consuming, and that it is 
sensitive to gradients in the disk density profile, and that 
using it would remove the simplicity of the J-grid approach 
we feel that it is a price worth paying. 

As we show below (in 331), and as discussed in Dutton 
et al. (2007), models with halo contraction (and standard 
stellar IMF's) are unable to reproduce the zero points of the 
VMR relations as well as the low galaxy formation efficiency 
required to reconcile the halo mass function and galaxy stel- 
lar mass function. While there are processes such as dynam- 
ical friction and impulsive mass loss that can expand the 
halo, implementing these in a galaxy evolution model is a 
non-trivial task. Thus for simplicity we wish to consider a 
model in which the dark halo does not respond to galaxy 
formation. 

Note that simply using the mapping between radius and 
j based on the halo profile (at each time step) in the absence 
of galaxy formation would not conserve the specific angular 
momentum of the baryons. To calculate the mapping be- 
tween j and radius, for the case of no adiabatic contraction, 
we solve the equation 



R = 



f/G 

M ha lo(< R) + M disk (< j) ' 



(19) 



where Mhaio(-R) is the mass (within a spherical radius R) of 
the dark matter plus hot gas halo in the absence of galaxy 
formation, and Mdisk(j) is the mass of the disk (gas plus 
stars) with specific angular momentum less than j. This 
way the self-gravity of the disk is included but adiabatic 
contraction is ignored. 



2.6 Star Formation 

Observations have shown that the disk averaged star forma- 
tion rates in nearby spiral galaxies are well fitted by a simple 
Schmidt (1959) law 



[M pc- 2 Gyr-i] [M pc~2 Gyr" 1 ] \[lM Q pc^] 



.(20) 



Kennicutt (1998) used S gas as the total gas density (molec- 
ular and atomic, but not including helium) and found esf = 
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10 20 30 1 2 3 

R [kpc] log 2 [M pc" 2 ] 

Figure 1. Right: Star formation surface density vs total gas surface density for our adopted star formation model [Eq. I|25jl ]. with 
egp = 13Gyr — . This has been calculated using the molecular and total gas densities of the model in the left panel. The dotted line 
shows the standard SK relation, which our model converges to at high gas density [Eq. H29H 1. The dashed line shows the relation our 
model converges to at low gas density Eq. H28I I. Left: Surface density profiles of stars and gas (total, atomic and molecular) which is used 
to calculate the right panel. The molecular fraction is computed using Eqs. 1221240 . 



0.25M Q p<r 2 Gyr"\ and n = 1.4 ±0.15 (Kennicutt 19980 
This simple empirical law holds over many orders of magni- 
tude in gas surface density, and even applies to circumnu- 
clear starburst regions. However, when applied to local gas 
densities, the Schmidt law breaks down at low gas densities, 
corresponding to large radii, where star formation has been 
found to be abruptly suppressed. Kennicutt (1989) argued 
that this suppression is due to the gas density falling below 
the critical density for global disk stability as given by the 
Toomre criterion (Toomre 1964) 



Scrit = a gas n(R)/3.36GQ, 



(21) 



where cr ga s is the velocity dispersion of the gas, n is the 
epicyclic frequency, and Q is the Toomre Q parameter. Ken- 
nicutt (1989) found that cr gas = 6kms _1 and Q = 1.5 repro- 
duces the observed star formation truncation radii. However 
other authors argue that this is just a coincidence, and that 
E cr it is not a threshold density (e.g. Schaye 2004). 

The physical origin of the Schmidt-Kennicutt (SK) re- 
lation is also not clear. However, something that is well es- 
tablished is that stars form out of molecular gas, and pre- 
dominantly in giant molecular clouds (GMC's). This lead 
Wong & Blitz (2002) to argue for a Schmidt law based on 
the surface density of molecular gas. For high gas densities 
the molecular gas dominates, so the two prescriptions are 
identical. However, for low gas densities the molecular frac- 
tion is suppressed, resulting in a steep dependence of the star 
formation rate density on total gas density. This gives an al- 
ternative explanation for the suppression of star formation 
at low gas densities. 



Including a helium correction of 1.36 results in egF 
0.16 M pc _2 Gyr -1 . 



The fraction of gas that is molecular, / mo i, can be de- 
fined in terms of the mass ratio between molecular and 
atomic gas, 7£ mo i by 



/moi 



Tin 



72-mol + 1 



(22) 



Empirically Wong & Blitz (2002) and Blitz & Rosolowsky 
(2004; 2006) have argued that 7£ mo i is determined to first 
order by the mid plane pressure, P ox t- The most recent re- 
lation from Blitz & Rosolowsky (2006) is 



Eaton 



f ext / k 



4.3 ±0.6 x 10 4 



(23) 



where k is Boltzmann's constant, and P e xt/fc is in cgs units. 
For a gas plus stellar disk the mid plane pressure is given, 
to within 10% by (Elmegreen 1993) 



— GE g 

2 E 



+ 



/ <7ga S \ 

V cr star / 



(24) 



where <T ga s and crstar are the velocity dispersions of the gas 
and stellar disk, respectively. For simplicity we will assume 
Cgas/o"star = 0.1 which is a reasonable assumption for the 
stellar dominated parts of disk galaxies. In the outer parts 
of disk galaxies this ratio is likely to be higher. However, 
in these regions gas dominates, and so the contribution of 
the stars to the mid plane pressure is negligible, regardless 
of the ratio between cr gas and <7 s t ar . In the inner regions of 
galaxies o" gaa /o" s t ar is likely to be smaller than 0.1, but here 
the densities are high enough that / m oi — 1- 

Following Blitz & Rosolowski (2006) we assume that 
star formation takes place in dense molecular gas, traced by 
HCN, with a constant star formation efficiency 



Esi 



[M pc- 2 Gyr- 



[Gyr-i] [M pc- 



(25) 
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where i S F ~ 10 - lSGyr" 1 (Gao & Solomon 2004, Wu 
et al. 2005). Expressing this equation in terms of the total 
gas content: 

Esfr ?sf E gi 



[M Q pc- 2 Gyr-i] [Gyr~i] [M pc" 2 ] 



/mol 72.HCN, 



(26) 



where 72-hcn = Emoi,HCN/E mo i is the ratio between the 
dense molecular gas (as traced by HCN) and the total molec- 
ular gas. 

Based on the arguments and references in Blitz & 
Rosolowski (2006) we adopt the following relation for 7?.hcn 



TZhcn = 0.1 1 



[200 M pc- 2 



(27) 



In the low pressure, low molecular density regime, 72-hcn 
0.1, and thus Eq. ([25} asymptotes to 



0.013 



[M Q pc- 2 Gyr-i] [Gyr-*] [M pc" 



[M pc- 



In the high pressure, high molecular density regime, 7?.hcn oc 
Emoi> and eq. (|25|l asymptotes to the familiar SK relation 



jSFR 



£sf 



0.012 



[M pc- 2 Gyr- 1 ] [Gyr-i] [M pc- 



[M pc- 



.(29) 



Furthermore, with ?sf = 13 Gyr , we recover the coefficient 
of esF = 0.16 of the standard SK relation. 

Fig. [T] shows the relation between star formation rate 
density and gas density for our star formation model 
Eq. (]25p - Note that in order to compute Esfr one needs to 
know E s tar and E gas . For illustrative purposes we have cho- 
sen a stellar and gas density profile representative of bright 
nearby disk galaxies (left panel of Fig. [T}. At small radii, 
and high gas densities / mo i — 1, and the star formation law 
follows the standard SK relation. At large radii the molec- 
ular fraction is very low, which results in a steeper slope to 
the star formation law. 

We implement the star formation recipe given by 
Eq. (|25[) as follows. At each time step and annulus in the 
disk, we calculate the star formation rate. Then we use the 
following approximation (valid for times steps small com- 
pared to the star formation time scale) to calculate the mass 
of newly formed stars 

AMstar(-R) = A(R) Esfr(-R, t) At, (30) 
with A the area of the annulus, and At the time step interval. 



2.7 Supernova Feedback 

When stars evolve they put energy back into the inter stel- 
lar medium. The effect of this on the star formation rate is 
partially taken into account by our empirically determined 
star formation recipe. What is not taken into account is a 
feedback driven outflow of gas from the disk. The physi- 
cal mechanism responsible for driving outflows is a subject 
of debate (e.g. Finlator & Dave 2008), so in this paper we 
consider both energy and momentum driven winds. 

Following van den Bosch (2001) we assume that the 
outflow moves at the local escape velocity of the disk-halo 
system. This choice is motivated by the fact that it maxi- 
mizes the mass loss from the disk-halo system (lower velocity 



winds will not escape the halo, and higher velocity winds will 
carry less mass). 

For our energy driven wind model following Dekel & 
Silk (1986) we assume that the kinetic energy of the wind is 
equal to a fraction, eefb, of the kinetic energy produced by 
SN. However, contrary to Dekel & Silk (1986) we apply this 
energy condition locally in the disk as a function of radius, 
rather than globally to the whole galaxy. Thus the mass 
ejected from radius, R, during a given time step is given by 



AM eject (R) = 



2 eefb Esn »?sn 
VL(R) 



(31) 



Here AM stal (iJ) is the mass in stars formed at radius, R, 
E S n = 10 51 erg ~ 5.0 x 10 7 km 2 s~ 2 M is the energy pro- 
duced by one SN, and ??sn = 8.3 x 10 -3 is the number of SN 
per solar mass of stars formed (for a Chabrier IMF). 
The local escape velocity is given by 



.(28) V CSC (R) = v/2|$ tot (i?)| 



(32) 



where <f>tot is the sum of the potentials due to the disk (stars 
plus gas) and halo (dark matter plus hot gas), and is com- 
puted assuming spherical symmetry. 

For our momentum driven wind model we assume that 
the momentum of the wind is equal to a fraction, cmfb, of 
the momentum produced by SN, thus the mass ejected from 
radius, R, during a given time step is given by 



A M 



fMFB PSN VSN 

V esc (R) 



AMstar(-R). 



(33) 



Here psn = 3 x 1O 4 M kms -1 is the momentum produced 
by one SN, assuming that each SN produces ~ 1OM of 
material moving at v ~ 3000 kms -1 (Murray, Quataert & 
Thompson 2005). Note that this corresponds to a kinetic 
energy of 4.5 x 10 7 M km 2 s~ 2 ~ 10 51 erg. 

We assume that the ejected mass is lost forever from the 
system: the ejected mass is not considered for later infall, 
and the corresponding metals are not used to enrich the 
infalling gas. This is clearly a dramatic oversimplification, 
but we make this choice to maximize the amount of gas 
that is lost from the system. 



2.8 Stellar Populations and Chemical Enrichment 

In order to convert stellar mass into luminosities we use 
the Bruzual & Chariot (2003) stellar population synthesis 
models. These models provide the luminosities L(t, Z) of a 
single burst stellar population with a total mass of 1M 
as a function of age, t, and metallicity, Z, in various optical 
pass-bands. To compute the luminosities of our model stellar 
populations we convolve the star formation history of our 
galaxies with the single burst stellar population synthesis 
models. 

In order to model the chemical enrichment of the ISM, 
we adopt the instantaneous recycling approximation. We as- 
sume that a fraction 1Z of the mass in stars formed is in- 
stantaneously returned to the cold gas phase with a yield y 
(defined as the fraction of mass converted into stars that is 
returned to the ISM in the form of newly processed metals) . 

The equations for the change in the cold gas mass and 
metals are: 



AM a 



AMcooi - (1 - K)AM slar - AM e j 



eject 



(34) 
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AAf mcta l — Zhot AiWcool — ^cold (1 - K)AM a 

+2/AM star - Z co i d AMejeot- (35) 

The metallicity of the cold gas is then given by Z co id = 
-Mmetai/AZcoid- Note we assume that the ejected metals do 
not enrich the hot halo gas. 

2.9 Book Keeping 

We now briefly describe how we keep track of the evolution 
of the various mass components. Given a, z — halo mass 
(dark matter plus hot gas) and concentration we compute 
the MAH of the halo using Eq. flTJ). The evolution of the 
virial radius and internal structure of the halo is then deter- 
mined by Eqs. ([3} & (fT0|) . We set up a grid in radius from 
0.001 to 1 times the redshift zero virial radius. As described 
in § 12.51 for the purposes of conserving angular momentum 
it is more convenient to use a grid in specific angular mo- 
mentum, j. Thus we convert the grid in R, to a grid in j 
using f/G = RM(R). 

For each time step (t — At, t) we compute the halo mass 
that is added to each radial bin so that the total mass follows 
the NFW profile for a halo of a given c(z) and Af vlr (2). Thus 

AM vir (j, t) = M vil (j, t) - Af vir (j, t-At). (36) 

We assume that the baryons make up a fraction /b ar of this 
mass, so that 

AMdmCj, *) = (!- /b ar )AM vir (j, t), (37) 

AM hoi (j, t) = / bar AA/virO', t). (38) 

When we compute the circular velocity we assume that the 
hot gas follows the mass distribution of the dark matter. 
When computing the cooling time we assume phot is /bar 
times the density at the virial radius at the time when the 
gas virialized. 

At each time step we then compute, using the recipes 
in the previous sections, the amount of gas that cools, 
AM coo i(j,i), the amount of stars formed AAf sta rO, i), and 
the amount of ejected gas AM c j oct (j,t). For the stellar pop- 
ulation modeling we keep track of the mass of stars formed 
at each time step and the metallicity of the gas from which 
the stars formed, Z co id(j, t). 

2.10 Overview of Model Parameters 

The input parameters of our models are as follows. 

(1) Cosmology: f2 m ,o, Ha, fib, erg, h, n. In this paper 
we adopt a flat ACDM cosmology motivated by the 5th 
year WMAP results (Dunkley et al. 2009), with fi m , = 
0.258, fi A = 0.742, fi b = 0.044, h = 0.7, cr 8 = 0.80 and n = 1. 

(2) Halo structure: K,F,ai nc . We adopt the Bullock 
et al. (2001a) model with F = 0.01, K = 3.7, and a\ nc = 
0.25. These parameters reproduce the distribution of halo 
concentrations of relaxed dark matter haloes in numerical 
simulations (Wechsler et al. 2002; Maccio et al. 2007, 2008). 

(3) Angular momentum distribution: A, /J\ n \, ot, oi oga . 
As fiducial values we adopt a median spin parameter A = 
0.035 with a scatter <j\ n \ = 0.54, corresponding to relaxed 
haloes (Maccio et al. 2007, 2008). For the angular momen- 
tum shape parameter we adopt a median a = 0.90, and 
scatter erioga = 0.11 (Sharma & Steinmetz 2005). 



(4) Star formation: esf- We use a star formation model 
based on dense molecular gas [Eq. (|25[) ]. and adopt esf = 
13 Gyr -1 . 

(5) Feedback: e E FB, £mfb, »?sn, JSsn, Psn- We adopt 
E SN = 10 51 erg, psn = 3xl0 4 Af kms" 1 , r/ SN = 8.3xl0" 3 . 
We treat ebfb and cmfb as free parameters. 

(6) Stellar populations and chemical enrichment: 
1Z, y, Zhot, and the choice of initial mass function (IMF). We 
adopt the Chabrier (2003) IMF, a return fraction K = 0.35, 
a stellar yield y = 0.02, and a metallicity of the hot gas of 
Z hot = 0.002 (~ 0.LZ o ). 



2.11 Overview of Output Parameters 

The output parameters of our models, that we discuss in 
this paper, are as follows: 

• Af v i r , total mass inside virial radius [Mq] 

• Vvir, circular velocity at the virial radius [kms -1 ] 

• Vmax.h, maximum circular velocity of the halo without 
galaxy formation [kms -1 ] 

• A/g a i, galaxy mass (stars and cold gas) [Mq] 

• A/ S tar, stellar mass [Mq]. 

• A/coid, cold gas mass [Af©]. 

• Zcoid, metallicity of cold gas. 

• €qp = m ga i/(fib/fim,o), g a l ax y formation efficiency. 

• figai = Afgai/Afvir, galaxy mass fraction. 

• fistar = Af s t a r /Af v i r , stellar mass fraction. 

• cricoid = Af co id /Mvir , cold gas mass fraction. 

• /gas = Afcoid/ (Afcoid + A/ S tar), cold gas fraction. 

• V2.21, the circular velocity measured at 2.15 /-band disk 
scale lengths [kms -1 ]. 

• 7?di, /-band disk scale length [kpc]. 

• jgai = Jgai/Jvir, galaxy angular momentum fraction. 

• jstar = Jstar/Jvir, stellar angular momentum fraction. 

• jcoid = Jcoid/Jvir, cold gas angular momentum fraction. 

• A g ai = A(j g ai/m ga i), galaxy spin parameter. 

• A SI ar = A(j s tar/m st ar), stellar spin parameter. 

• A co id = A(j co id/m co id), cold gas spin parameter. 

• A/bar = /bar Afvir, the baryonic mass accretion rate 
[M Q Gyr -1 ] 

• A/ coo i, the mass cooling rate [Mq Gyr -1 ] 

• Msf, the star formation rate [Mq Gyr -1 ] 

• Af w , the mass outflow rate [Mq Gyr -1 ] 

• rj — M-w/Mst, the mass loading factor 

• Vw = £iV w (i)A/ w (i)/Af w , the (outflow) mass weighted 
mean outflow velocity [kms -1 ], where the i refers to the 
radial grid position. 

The disk scale lengths are determined using the follow- 
ing procedure, which was developed to give robust disk scale 
lengths for the full range of surface brightness profiles pro- 
duced by our model. We first compute the local disk scale 
length between the radii enclosing 50 and 99% of the stel- 
lar mass. The local disk scale length is computed at radial 
bin i by using the surface densities and radii at radial bins 
i — 1 and i + 1. We then determine the maximum of the 
local disk scale length, and the radius where this maximum 
occurs, -Rmax- Finally we determine the scale length of the 
disk using a linear fit to the model data over the range 
0.67? max < R s£ 1.67? max ■ 
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2.12 Comparison with other Disk Galaxy 
Structure Models 

In this section we place our model in the context of exist- 
ing analytic and semi-analytic models of disk galaxy struc- 
ture/formation in the literature. We classify these models 
into two general types: 1) models that conserve total spe- 
cific angular momentum (i.e. the structural profile of the 
disk is assumed) and 2) models that conserve the distribu- 
tion of specific angular momentum (i.e. the structural profile 
of the disk is derived). Both of these classes of models can 
be static or include evolution. The essential assumption in 
both classes of models is that the disk is in centrifugal equi- 
librium inside some potential (which may or may not include 
the self gravity of the disk). 



2.12.1 Models that Conserve Total Specific Angular 
Momentum 

In the simplest models of this class, the circular velocity is 
assumed to be constant, i.e. corresponding to an isothermal 
density profile, and the self gravity of the disk is ignored. 
This model has 3 parameters: the circular velocity, the spin 
parameter and the disk mass fraction. Such a model was 
discussed in Mo, Mao, & White (1998; MMW), and is widely 
used in Semi Analytic Models (e.g. Kauffmann, White, & 
Guiderdoni 1993; Cole et al. 1994; Somerville & Primack 
1999; Hatton et al. 2003; Croton et al. 2006). 

A more realistic version of this model includes the self 
gravity of the disk and adiabatic contraction of the halo 
(Blumenthal et al. 1986), which usually results in smaller 
sizes for a given spin parameter and disk mass fraction. In 
Mo, Mao & White (1998), the halo was assumed to be an 
NFW profile, and the disk was assumed to be exponential. 
This model has 4 parameters: the circular velocity of the 
halo, the concentration of the halo, the spin parameter and 
the disk mass fraction. This version of the MMW model is 
widely used in Semi Analytic Models (Cole et al. (2000); 
Benson et al. (2003); Somerville et al. 2008) and studies of 
disk galaxy scaling relations (e.g. Navarro & Steinmetz 2000; 
Pizagno et al. 2005; Dutton et al. 2007; Gnedin et al. 2007). 



2.12.2 Models that Conserve the Distribution of Specific 
Angular Momentum 

The MMW type model is useful for understanding the origin 
of disk galaxy scaling relations, but it does not explain the 
origin of the density profiles of galaxy disks, or the relation 
between gas and stars in galaxy disks. In order to address 
these questions, one needs to start from some specific an- 
gular momentum distribution (AMD). This AMD may be 
that of a sphere in solid body rotation, or preferably that 
found in cosmological simulations (e.g. Bullock et al. 2001b). 
The radial density profile of the disk is then determined by 
detailed conservation of specific angular momentum. 

As with the MMW type models, these models may ig- 
nore the self gravity of the disk by assuming the total den- 
sity profile is isothermal (e.g. Ferguson & Clarke 2001), or 
include the self gravity of he disk inside a dark matter halo 
(e.g. Dalcanton, Spergel & Summers 1997). These models 
may also include evolution of the baryonic disk, by follow- 
ing the cooling of gas inside growing dark mater haloes, and 



evolution of the stellar disk, by following the star formation 
as a function of radius (e.g. Firmani & Avila Reese 2000; 
van den Bosch 2001, 2002; Stringer & Benson 2007). 

Our model fall into this later category, being an evolu- 
tion of the van den Bosch (2001; 2002) models. The strength 
of these models over the MMW type models is that they can 
be used to self-consistently study the origin and evolution 
of disk density profiles (stars, gas, star formation rates, in- 
flows, outflows, metallicity, stellar ages) and rotation curve 
shapes. 



3 THE DEPENDENCE OF THE GALAXY 

VELOCITY - MASS - SIZE RELATIONS ON 
FEEDBACK 

In this section we investigate the dependence of the scaling 
relations between rotation velocity, stellar mass and stellar 
disk scale size on the feedback model. 



3.1 Observed Disk Galaxy Scaling Relations 

Here we overview the main observed velocity, mass, size, 
scaling relations that we are going to compare our models to. 
We use the relations between rotation velocity at 2.2 /-band 
disk scale lengths, V2.21, stellar mass, M star , and /-band disk 
scale length, Rdi, from the data set of Courteau et al. (2007), 
as presented in Dutton et al. (2007). The stellar masses in 
Dutton et al. (2007) were derived from /-band luminosities 
using the relations from Bell et al. (2003) corresponding to 
a diet Salpeter IMF. Here we adopt a Chabrier IMF, and so 
subtract 0.10 dex from the stellar masses. 
The stellar mass TF relation is given by 

log , V2 - 2 \, = 2.195 + 0.259 (log , M ^ r - - 10.5 ) , (39) 

the size-stellar mass relation is given by 

[h^kpc] 



log . = 0.491 + 0.281 log - 10.5 ) , (40) 

[h 70 M Q \ 



and the corresponding size- velocity relation is given by 



log 



Rdi 



= 0.491 + 1.086 log 



V2.2 



[kms 



2.195 .(41) 



The intrinsic scatter in these relations is estimated to be 
Clog 10 v|Af — 0.05, o- logloJ j| A/ ~ 0.13, and (J logloii | V ~ 0.16. 
The errors on the slopes of the VM, RM, and RV relations 
from fitting uncertainties are 0.01, 0.02, 0.12, respectively. 
However, systematic uncertainties are significantly larger, 
and harder to quantify. The most significant selection effect 
for the slope of the RM relation is surface brightness. The 
data set compiled by Courteau et al. (2007) is likely missing 
lower surface brightness galaxies, and thus over-estimates 
the slope of the size-mass relation at low masses. Such a 
conclusion is supported by Shen et al. (2003) who studied 
the half light radius-stellar mass relation (-R50 — iWatar) for 
a much larger sample of galaxies (~ 10 ) from the SDSS. 
They find a log slope of 0.14 at low masses, increasing to 
0.39 at high masses. 



© 2009 RAS, MNRAS OOO.mWl 



10 Button & van den Bosch 




Figure 2. Effect of energy feedback efficiency, eEFBi on the position of a galaxy with M v i r = 6.3 X lO 11 ^ - 1 Mq in the VMR planes. 
Models have feedback efficiencies of £efb = 0,0.1,0.25,0.5, and 1.0. The arrows indicate the direction of increasing eefb- Models with 
adiabatic contraction are shown with solid red symbols, models without adiabatic contraction are shown as black open symbols. In the 
V M plane the horizontal dotted grey lines show virial velocity of the halo, V viv , and the maximum circular velocity of the halo prior to 
galaxy formation, V maXiri . The solid and dashed green lines in the VMR panels show the mean and 2a scatter of the observed relations 
from Dutton et al. 2007 (D07), assuming a Chabrier IMF. The long dashed green line shows the observed half-light radius stellar mass 
relation from Shcn et al. (2003). The panel in the top right shows the effect of feedback on the galaxy mass fraction, m ga i (circles), 
and galaxy spin parameter, A ga i (triangles). The dashed horizontal lines show galaxy formation efficiencies of 100, 50, and 25 percent, 
the dotted horizontal line shows the spin parameter of the halo. As the feedback efficiency is increased the galaxy mass fraction (m ga i) 
decreases, the galaxy spin parameter (A ga i) increases, the rotation velocity decreases, the stellar mass decreases, and the size of the stellar 
disk increases. 



3.2 A Fiducial Model 

To illustrate the effect that feedback has on observable prop- 
erties of disk galaxies we consider a model with virial mass 
Mvir = 6.3 x 10 n /i -1 M©, and median concentration and 
angular momentum parameters: c v i r = 9.9, A = 0.035 and 
a = 0.9. The results of varying the energy and momentum 
feedback efficiency parameters from to 1 are shown in Figs. 
[5] & El respectively. The main panels show the V M R rela- 
tions with the solid line showing the mean relations from 
observations and the dashed lines showing the 2<r scatter. 
The upper right panels in these figures shows the more the- 
oretical parameters m ga i and A ga i. These are not directly ob- 
servable because they require knowledge of the halo mass, 
a quantity that cannot, at present, be reliably measured for 
individual galaxies. 

We first focus on the models with energy feedback and 
adiabatic contraction (solid red points and lines in Fig. [2}. 



A model with no feedback (eefb = 0) results in a galaxy 
mass fraction of ~ 0.8/bar (where /bar — 0.17 is the univer- 
sal baryon fraction). The galaxy formation efficiency is less 
than 1 because cooling starts to become inefficient at late 
times. Since the highest specific angular momentum mate- 
rial is accreted last, and this does not have time to cool, the 
spin of the galaxy is slightly lower than that of the halo. 
The high galaxy formation efficiency results in a disk scale 
length a factor of ~ 1.8 too small and a circular velocity at 
2.2 disk scale lengths a factor of ~ 1.5 too high. 

When feedback is included, some of the cold baryons 
are ejected from the disc (and halo). This trivially results 
in lower m ga i (hence lower baryonic mass) but also, non- 
trivially, higher A ga i. The higher A ga i is due to the prefer- 
ential ejection of low angular momentum material, which 
we discuss in more detail in fj4] Both of these effects re- 
sults in larger, lower surface density disks, which result in 
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Figure 3. As fig [2] but for momentum driven winds. 



less efficient star formation, and hence higher gas fractions 
and lower stellar masses. The reduction in the amount and 
density of the baryons results in a lower rotation velocity, 
both because the baryons contribute less to the circular ve- 
locity at 2.2 disk scale lengths and because there is less 
halo contraction. With a high energy feedback efficiency of 
cefb = 1-0 the galaxy formation efficiency drops to ~ 15%: 
~ 20% of the baryons have not cooled, while ~ 65% of 
the baryons have been ejected from the disk and halo. This 
model also has a galaxy spin a factor of ~ 1.7 higher than the 
halo spin. The low galaxy mass fraction and higher galaxy 
spin parameter results in sizes that are more than a factor of 
~ 2 too large. The feedback efficiency can be tuned so that 
the model reproduces the size-mass, or size-velocity relation. 
However, for all feedback efficiencies, the models rotate too 
fast at a given stellar mass. This is because higher feedback 
efficiency results in lower stellar masses as well as lower rota- 
tion velocities, with the net result that galaxies move almost 
parallel to the TF relation. 

Fig. [3] shows the same results as Fig. [2] but for the 
momentum driven wind model. Even with 100 per cent ef- 
ficiency this galaxy formation efficiency is ~ 50%, where 
~ 30% of the baryons have been ejected from the disk and 
halo. This results in sizes that are in agreement with obser- 
vations, but the models still rotate too fast. The reason that 
energy driven winds are more efficient at ejecting mass than 



momentum driven winds is shown in Fig. [4] which shows 
the mass loading factor, 77, versus the wind velocity, V w . 
The mass loading factor is defined as the ratio between the 
mass outflow rate and the star formation rate. For energy 
driven winds r\ oc V^ 2 , whereas for momentum driven winds 
77 oc V^T 1 . Thus everything being equal, energy driven winds 
have higher mass loading factors than momentum driven 
winds for all typical wind velocities 

3.3 The Tully-Fisher Zero Point Problem 

A common problem to both feedback models, for all values of 
the feedback efficiency, is that they over predict the rotation 
velocities. This is a standard problem for galaxy formation 
models in ACDM. As discussed in Dutton et al. (2007) and 
Gnedin et al. (2007), there are 3 solutions: 1) Lower the 
stellar mass-to-light ratio (i.e. for a given luminosity there is 
less stellar mass which shifts observed galaxies to the left and 
hence higher velocity in the VM plane jj. 2) Lower the initial 
halo concentration (which directly lowers V2.2). 3) Reverse 
halo contraction (which directly lowers V2.2). 

2 Note that the stellar masses of the model galaxies are not com- 
pletely independent of the IMF, as the return fraction is IMF 
dependent. However, changes in the return fraction are compen- 
sated for by changes in the star formation rate. 
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Figure 4. Mass loading factor versus wind velocity for energy 
and momentum driven winds. For equal efficiencies an energy 
driven wind is always more efficient than a momentum driven 
wind, especially for low wind velocities. 



A small change in the stellar mass-to-light ratio would 
be plausible due to systematic uncertainties (such as in the 
IMF or the stellar populations synthesis models) in the mea- 
surement of stellar mass-to-light ratios. However, the stellar 
mass-to-light ratios would have to be lowered by about a 
factor of 2 to match the TF zero point. Such a large change 
would require a top-heavy IMF. But all of the available con- 
straints on stellar mass-to-light ratios point to IMF's similar 
to Chabrier (e.g. de Jong & Bell 2007). 

Lower halo concentrations would require less power on 
galaxy scales than in standard ACDM. This would also re- 
duce the amount of substructure, which could help solve 
the missing satellite problem (Klypin et al. 1999; Moore 
et al. 1999). However, the recent discovery of many satel- 
lite galaxies around the Milky Way has lessened the dis- 
crepancy between observations and standard ACDM (e.g. 
Tollerud et al. 2008 and references therein). Furthermore us- 
ing cosmological simulations with parameters from the lat- 
est WMAP results (Spergel et al. 2007; Dunkley et al. 2009), 
Maccio et al. (2008) have shown that the central densities 
of dark matter haloes are consistent (both in normaliza- 
tion and scatter) with those measured from dwarf and LSB 
galaxies (which typically have maximum rotation velocities 
of ~ 100 km s -1 ). Thus there does not seem compelling evi- 
dence for a modification to ACDM on small scales. 

Given that reducing stellar mass-to-light ratios or the 
initial halo concentrations do not seem plausible, we consider 
the third possibility, that halos do not contract as expected. 
There are two processes that could cause the halo to expand. 
1) Dynamical friction between baryons and the halo, e.g. by 
infalling baryonic clumps (e.g. El Zant et al. 2001; El-Zant 
et al. 2004; Tonini et al. 2006; Jardel & Sellwood 2009); or by 
galactic bars (Weinberg & Katz 2002; Holley-Bockelmann, 
Weinberg & Katz 2005; Sellwood 2008); 2) Rapid (i.e. non- 



Table 1. Model Parameters 



Name 


e EFB 


£ MFB 


A 


0"lnA 


AC 


I: No Feedback: 


0.0 


0.0 


0.035 


0.35 


N 


II: Momentum: 


0.0 


1.0 


0.035 


0.35 


N 


III: Energy: 


0.25 


0.0 


0.025 


0.35 


N 



adiabatic) mass loss from the galaxy, e.g. by SN driven winds 
(Gnedin & Zhao 2002; Read & Gilmore 2005). Both of these 
effects have been shown to be effective at expanding the 
halo, but the combined effect (which may be greater than 
the sum of its parts) has so far not been investigated. Note 
that our standard model with adiabatic contraction takes 
into account the adiabatic expansion of the halo due to out- 
flows. However, since there is a net inflow of baryons into 
the centers of galaxies, the overall effect is halo contraction. 

Furthermore halo contraction is based on the idea that 
galaxies form by cooling flows. The hot gas radiates away its 
binding energy, so when it falls to the center of the potential 
to form the galaxy the halo has to contract. However, recent 
simulations have indicated that the gas, in haloes that host 
disk galaxies, is accreted by cold flows (i.e. it does not shock 
heat when it enters the virial radius) . This new scenario thus 
allows the possibility of the baryons transferring energy to 
the halo during galaxy formation. 

The black circles in Figs. [5] & [3] show the results for 
two models without halo contraction. For both energy and 
momentum driven winds the galaxy mass fractions and 
spin parameters have the same dependence on the feed- 
back efficiency in models with and without halo contraction. 
However, the models without halo contraction have signifi- 
cantly lower rotation velocities, allowing a match to the TF 
zero point for energy feedback efficiencies of ~ 0.1 — 0.5, 
or momentum feedback efficiencies of ~ 1. A model with 
£efb = 0.25 also has V2.21 — V maXi h, and a galaxy formation 
efficiency of ~ 35%, consistent with observations (see §1). 
However, this model has sizes that are too large. The solu- 
tion to this is to lower the spin parameter of the baryons. 
This could occur if disk galaxies formed in haloes with lower 
than average halo spin, or by the baryons transferring some 
of their angular momentum to the halo during galaxy as- 
sembly, for example via dynamical friction. As we discussed 
above, such a process may be responsible for expanding the 
halo. As we show below, a model with a factor 1.4 lower 
halo spin reproduces both VM and RM relations as well as 
low galaxy formation efficiencies. The momentum feedback 
model with maximum feedback efficiency, emfb = 1.0, on 
the other hand, roughly reproduces the VM and RM zero 
points, but it has galaxy mass fractions and V2.2i/Knax,h 
that are too high compared to observational and theoretical 
constraints. 



3.4 Models with Scatter 

Having discussed the effects of energy vs. momentum driven 
winds, and halo contraction vs no halo contraction for a 
single halo mass, with the median concentration, and an- 
gular momentum parameters, we now turn our attention to 
models with the full range of halo masses relevant for disk 
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Figure 5. V2.21 — A/star — Rdl relations for model I (no feedback) for virial masses logarithmically sampled in the range 10 10 < Af v i r < 
10 13 /i — 1 M.q. The grey points and black lines show the models while the green lines show the observed relations. In the VMR panels the 
solid and dashed green lines show the mean and la scatter, assuming a Chabrier IMF, of the observed relations from Dutton et al. 2007 
(D07). In the RM panel the long dashed green line shows the observed half-light radius stellar mass relation from Shen et al. (2003); In 
the VM — RM residual panel the green dashed line shows the observed correlation from Avila- Reese et al. (2008). The black lines show 
the 16th, and 84th percentiles of the models, while the error bars show the Poisson error on the median. The thick black line shows a 
double power law fit to the median using the function in Eq. 1142 II . the parameters of the fits are given in Table [2] The solid red lines show 
the mean of the models fitted over a range corresponding to 10 9 ' 6 < Af s t ar < 10 110 Mq . The parameters of the fits (zero point, slope, 
scatter (in log 10 units)) are given in the top left corner of each panel. The relations are fitted as follows: log V = a + 6 (log M — 10.5); 
log R = a + b (log M — 10.5); log R = a + b (log V — 2.195). The upper right panel shows the residuals of the VM relation vs the residuals of 
the RM relation. The red lines show the mean and la scatter of a fit of the form: A log V\M = b Alogi?|Af. The correlation coefficient, 
r, is also given where r = ba x /a y , where x = Alogi?|M and y = A log V\M. This model fails to reproduce all of the observations (with 
the exception of the slope of the VM relation). 



galaxies, and with distributions of concentration, and angu- 
lar momentum parameters. 

We run a Monte Carlo simulations with halo masses 
ranging from 10 10 < Af v i r < 10 13 h~ 1 M ( r ) , corresponding 
to virial velocities ranging from 31 < Vvir J$ 310 km s -1 . 
In ACDM there are many more low mass haloes than high 
mass haloes, however, since we are interested in the scaling 
relations between galaxies, rather than the number densities 
we sample the halo masses uniformly in log-space. 

As discussed in Dutton et al. (2007) we also find that 
models with the expected scatter in halo spin parameter 
finA = 0.54 significantly over predict the amount of scat- 
ter in the RM and RV relations. This may signify that 
disk galaxies form in a special sub-set of haloes, or that the 
baryons acquire a different distribution of specific angular 



momentum than the dark matter. For the remainder of this 
paper we adopt cr\ n x ~ 0.35, as this provides a reasonable 
agreement to the observed scatter in disk sizes. 

To illustrate the effect of feedback on galaxy scaling re- 
lations we consider three models. Model I has no feedback, 
Model II has momentum feedback, and Model III has en- 
ergy feedback and an average spin parameter a factor of 1.4 
lower than models I and II. The parameters of these models 
are given in Table [1] The parameters of Model II and III 
were chosen to match the zero points of the VMR relations, 
and thus for reasons discussed above, they have no adiabatic 
contraction. For each model we run a Monte-Carlo simula- 
tion consisting of 2000 galaxies. For each galaxy we select 
the parameters c, A, and a from log-normal distributions 



© 2009 RAS, MNRAS 000.1111271 



14 Button & van den Bosch 




Figure 6. As Fig. [5]but for Model II (momentum driven feedback). This model provides a reasonable match to the zero points, slopes 
and scatters of the V MR relations. 



with means and scatters as determined by the parameters 
in § l2~TTJl and Tabled] 

Figs. !5l7l show the VMR relations, as well as the corre- 
lation between the residuals of the VM and RM relations, 
for models I-III. Recall that V is the circular velocity mea- 
sured at 2.15 /-band disk scale lengths, M is the stellar 
mass, and R is the /-band disk scale length. 



3.4-1 slopes 

The VMR relations in these figures are fit with two relations: 
a single power-law over the range where there is observa- 
tional data, and a double power-law over the full range of 
masses. The parameters of the best-fit single power-law fits 
are indicated in the panels of Figs. 15171 The double power- 
law is given by 





T 1 




)■ 




(-)l 




.2 + 2 


\XoJ J 



Here a is the slope at x <C xq\ (3 is the slope at x ^> xq; xq 
is the transition scale; j/o is the value of y at xq. The best 
fit values of these parameters are given in Tabled 

The slope (as given by the single power-law fits) of 
the VM relation is only weakly dependent on the feedback 
model. This is expected since (as shown in Figs. [2] & [3} the 



Table 2. Parameters of double power-law fits to VMR 
relations in Figs. 15171 using the function in Eq. 1 142 I I . 



Model 


a. 





A/ a tar, 


V2.2/,0 


I 


0.27 


0.34 


10.71 


2.34 


11 


0.29 


0.44 


10.94 


2.36 


111 


0.22 


0.46 


10.35 


2.17 


Model 


a 


P 


A/star, 


R<a,a 


I 


0.44 


0.14 


11.56 


0.65 


11 


0.28 


0.17 


11.15 


0.68 


111 


0.19 


0.04 


10.55 


0.52 


Model 


a. 


P 


^2.21,0 


Rdl,0 


I 


1.79 


0.40 


2.59 


0.65 


11 


1.15 


0.00 


2.50 


0.72 


111 


1.15 


0.00 


2.00 


0.40 



offset of a galaxy from the VM relation is only weakly de- 
pendent on the galaxy mass fraction, which is determined 
by feedback (for haloes with M V \ T < 1O 12 M0). By contrast 
the slope of the RM relation depends strongly on the feed- 
back model. Again, this is expected given that the offset of a 
galaxy from the RM relation depends strongly on the galaxy 
mass fraction. The model without feedback (model I) has a 
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Figure 7. As Fig. [5] but for Model III (energy driven feedback). This model provides a reasonable match to the zero point, slope and 
scatter of the VM relation. It predicts a shallower slope to the RM and RV relations than our fiducial data from Dutton et al. 2007, 
however at low masses the shallow slope of the model is in good agreement with the data from Shen et al. (2003) (green long dashed). 
The correlation between the residuals of the VM and RM relations (top right panel) of the models is stronger than the observations, 
although observational uncertainties in stellar mass measurements will cause the observed correlation to be underestimated. 



slope of 0.40, the model with momentum driven feedback 
(model II) has a slope of 0.26, and the model with energy 
driven feedback (model III) has a slope of 0.14. 

The observed slope of the size-mass relation from 
(Courteau et al. 2007 and Dutton et al. 2007) is ~ 0.28, 
which favors the momentum driven wind model. However, as 
discussed in § 13.11 at low stellar masses (M star < 10 10 M Q ), 
this slope is likely biased high by selection effects. Shen 
et al. (2003) find a slope of 0.14 at low masses for the 
half-light radius-stellar mass relation for a much more com- 
plete sample of late-type galaxies. Such a shallow slope is 
in much better agreement with our energy driven model. 
Thus the observed slope of the RM relation favors the en- 
ergy driven wind model at low masses, and the momen- 
tum driven model at high masses. However, at high stel- 
lar masses (M sta r ^ 1O 1O ' 5 M0) bulges are common is spiral 
galaxies (e.g. Weinzirl et al. 2008). Due to the correlation 
between the masses of bulges and black holes (Magorrian 
et al. 1998), AGN feedback may play a significant role in 
regulating galaxy formation efficiency in high mass spiral 
galaxies. Thus it is plausible that a model with energy driven 
SN feedback (which primarily effects galaxies in low mass 



haloes) and AGN feedback (which primarily effects galaxies 
in high mass haloes) could explain the slopes of the size- 
mass relation at low and high masses. However, since AGN 
feedback is not expected to be significant for galaxies in low 
mass haloes, it is unlikely that AGN feedback will be able 
to help the momentum driven wind model reproduce the 
shallow slope of the size-mass relation at low masses. 



3.4-2 scatter and residual correlations 

All models produce a VM relation with relatively small scat- 
ter, with smaller scatter in the models with feedback. The 
amount of scatter in the VM relation is directly related to 
the strength of the correlation between the VM and RM re- 
lations. The model without feedback has a very strong corre- 
lation (correlation coefficient, r = —0.98; slope, b = —0.39), 
which is caused by these galaxies being baryon dominated 
at 2.2 disk scale lengths. The models with feedback have 
weaker, but still significant, correlations (r ~ —0.8, b ~ 
—0.2). These correlations are stronger than that observed 
by Courteau et al. (2007) for the /-band VL and RL rela- 
tions (r = -0.16, b = -0.07), and by Gnedin et al. (2007) 
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Figure 8. Dependence of the Baryonic Tully-Fisher (BTF) relation at rcdshift z = on feedback. The green points show data from 
McGaugh (2005), where the rotation velocities have been measured in the "flat" part of resolved HI rotation curves. Green lines show 
the mean and ltr scatter of a fit to this data. The grey points show the models, where VgOc is the circular velocity at a radius enclosing 
80% of the gas. The black lines show the 16th and 84th percentiles in bins of width 0.25 dex in M ga i. The red lines show a fit of the 
form: log Vgo c = a + 6(logAf ga i — 10). The fits are performed over the range that the red lines are plotted. The zero point, slope, and 
scatter for each fit are given in the upper right corner of each panel. The slope, zero point and scatter of the BTF relation is only weakly 
dependent on the feedback model. 



(r = 0.23 ± 0.14) and Avila- Reese et al. (2008) (r = -0.29, 
b — —0.09) for the stellar mass VM and KM relations. 
However, as discussed in Dutton et al. (2007) scatter in the 
stellar mass-to-light ratio of ~ 0.15 dex, either from intrin- 
sic variations or measurement uncertainties will weaken the 
correlation between the observed relations compared to the 
theoretical VM and RM relations. Thus we do not consider 
this failure of the model as a serious shortcoming. 



3.5 The Baryonic Tully-Fisher Relation 

The fundamental basis of the Tully-Fisher (linewidth- 
luminosity) relation is believed to be the relation between 
the asymptotic rotation velocity of a galaxy disk, Va at , and 
the baryonic mass, M ga i, (the sum of stellar and cold gas 
mass). This relation is referred to as the Baryonic Tully- 
Fisher (BTF) relation. It was first studied by McGaugh 
et al. (2000), and subsequently by Bell & de Jong (2001), 
McGaugh (2005), Geha et al. (2006), Noordermeer & Ver- 
heijen (2007), and Avila- Reese et al. (2008). 

The most significant source of uncertainty in the BTF is 
how one measures stellar mass. McGaugh (2005) measured 
the BTF for stellar masses calculated using different meth- 
ods: stellar population synthesis models (e.g. Bell & de Jong 
2001), the maximum disk hypothesis (van Albada & San- 
cisi 1986), the mass-discrepancy acceleration relation (i.e. 
Modified Newtonian Dynamics, or MOND, Milgrom 1983)). 
McGaugh (2005) found that the scatter in the BTF was 
minimized when the stellar masses were calculated with the 
mass-discrepancy acceleration relation. Under the assump- 
tion that the correct method would minimized the scatter 
in the BTF, this is evidence in favor of MOND over dark 
matter. 

However, this is a circular argument because a rela- 
tion between the asymptotic rotation velocity of a galactic 
disk, Vfl at and the baryonic mass, M ga i, with zero scatter, 
is built into MOND. Thus if the stellar masses are chosen 
based on the MOND prescription, they will result in a BTF 



with scatter only due to measurement errors on Vfl a t and 
distance uncertainties. The scatter in the BTF (as defined 
as the relation between Vk a t and M ga i) thus cannot be used 
to discriminate between MOND and dark matter. However, 
MOND generally predicts higher stellar masses than stel- 
lar population synthesis models (based on a Kroupa IMF). 
Thus if stellar masses could be measured independently, this 
would provide a means of falsifying MOND. 

Here we use the data from McGaugh (2005), using the 
stellar population synthesis stellar mass-to-light ratios, with 
an offset of -0.1 dex (corresponding to a Chabrier IMF). 
The majority of galaxies in this sample are in the UMa Clus- 
ter, for which the distance is somewhat uncertain. McGaugh 
(2005) adopted a distance of 15 Mpc. We adopt the HST Key 
Project distance of D = 20.7 Mpc (Sakai et al. 2000), which 
is also the distance used by Bell & de Jong (2001). 

The BTF data are plotted as green filled circles in Fig. [8] 
A linear fit to the data gives the following relation between 
the rotation velocity and baryonic mass 



log- 



kmt 



2.027 + 0.279 log 



gal 



10 



(43) 



with a scatter of 0.053 dex in logVfl at . This BTF is con- 
sistent with that from Bell & de Jong (2001) who report a 
slope of 0.285(±0.015) and a zero point of 2.031(±0.011), 
and that from Avila-Reese et al. (2008) who report a slope 
of 0.306(±0.012) a zero point of 2.036(±0.129) and an in- 
trinsic scatter of 0.051 dex in log V. The slightly steeper 
slope obtained by Avila-Reese et al. (2008) can be attributed 
(see Verheijen 2001) to their use of Hi linewidths, com- 
pared to Vfl a t used by McGaugh (2005) and Bell & de Jong 
(2001). This good agreement is reassuring given that the 
data samples are largely based on the data set of Verheijen 
et al. (2001). 

The BTF relations for our models are given by the grey 
dots in Fig. [8] For the rotation velocity we use Vsoc, the 
circular velocity at a radius enclosing 80% of the gas mass, 
Rsoc, which usually corresponds to the flat part of the ro- 
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tation curve (see Dutton 2008 rotation curves). Power-law 
fits to the models over baryonic masses between 1 x 10 9 and 
1 x 10 11 are shown as red lines in the figure. The param- 
eters of these fits are given in the top left of each panel. 
All three models result in BTF relations with similar slopes, 
zero points and scatter, and in reasonable agreement with 
observations. We note that for the model galaxies the slope 
of the BTF depends on galaxy mass, with slightly larger 
slopes for higher mass galaxies. 

The BTF relation has been used to constrain the rela- 
tion between baryonic mass and halo mass. By comparing 
the observed slope of the BTF (0.27 ± 0.01) to the predic- 
tion from CDM (the slope of the V ma:x ,h — M V [ T relation for 
dark matter haloes is 0.294 ± 0.005, Bullock et al. 2001a), 
Geha et al. (2006) argued that low mass galaxies have not 
preferentially lost baryons as would be predicted by feed- 
back models (e.g. Dekel & Silk 1986). However, this is based 
on the incorrect assumption that the maximum observed 
rotation velocity is equal (or proportional) to the maxi- 
mum circular velocity of the halo, independent of the baryon 
to halo mass ratio. As discussed by several authors (e.g. 
Navarro & Steinmetz 2000; Dutton et al. 2007; Avila-Reese 
et al. 2008), the maximum rotation velocity of a galaxy is, 
in general, not equal to the maximum circular velocity of 
the halo in the absence of galaxy formation, Knax,h- As the 
baryon fraction increases, so to does the maximum circu- 
lar velocity. This is because the baryons contribute a non- 
negligible amount of mass to the circular velocity. Thus for 
reasonable galaxy mass fractions, variation in galaxy mass 
fraction moves galaxies roughly parallel to the BTF. 

Fig. [8] shows that the slope, zero point and scatter of 
the BTF are remarkably insensitive to the feedback model. 
Furthermore, as shown in §4.1, our energy feedback model 
results in substantial differential mass loss between haloes 
of mass M vir ~ 10 10 M Q and M vir ~ 1O 12 M . Yet it has 
the same BTF slope as a model with no mass loss (and 
constant baryon to dark matter ratio within this range of 
halo masses). This provides a counter example to the claim 
by Geha et al. (2006) that models with preferential mass 
loss in dwarf galaxies cannot explain the slope of the BTF. 



4 THE DEPENDENCE OF GALAXY MASS 
FRACTIONS AND SPIN PARAMETERS ON 
FEEDBACK 

The two main parameters that determine the structure of 
disk galaxies are the mass and specific angular momentum 
of the cold gas and stars. We now investigate the dependence 
of these parameters on feedback and halo mass. 

4.1 Galaxy Mass Fractions 

Fig. [9] shows the mass fractions and spin parameters of our 
3 models as a function of virial mass. 

The relations in Fig.|5]are fitted with the following dou- 
ble power-law: 

/ x \ a Tl 1 / x \7-i(P-«)/7 
» = *U) L2 + 2UJ \ (M) 
Here a is the slope at x -C xo; /3 is the slope at x S> xo; 
xq is the transition scale; yo is the value of y at xq; and 7 



Table 3. Parameters of double power-law fits to mass 
fractions and spin parameters vs virial mass in Fig. [9] 
using the function in Eq. 11441 . 



Model 


a 


P 


7 


Mvir 


m cold,0 


I 


-0.06 


-0.57 


2.00 


11.85 


-1.58 


11 


-0.07 


-0.57 


2.45 


11.84 


-1.58 


111 


0.05 


-0.59 


1.50 


11.80 


-1.73 


Model 


a 


/3 


7 


-^vir.O 


m s tar,0 


I 


0.02 


-0.45 


3.54 


11.96 


-0.93 


11 


0.25 


-0.33 


3.60 


11.92 


-1.19 


111 


0.65 


-0.27 


2.55 


12.09 


-1.32 


Model 


a 


P 


7 






I 


0.00 


-0.51 


2.59 


12.02 


-0.86 


11 


0.10 


-0.43 


2.71 


12.01 


-1.06 


111 


0.26 


-0.29 


3.90 


12.14 


-1.19 


Model 


a 


P 


7 


A^ v j r 


A GO i(j 


I 


0.01 


-0.35 


4.14 


12.02 


-1.14 


11 


0.03 


-0.37 


2.45 


12.03 


-1.15 


111 


0.01 


-0.37 


2.95 


12.03 


-1.26 


Model 


a 


P 


7 




■^star 


I 


0.04 


-0.34 


2.79 


11.97 


-1.72 


11 


0.01 


-0.33 


3.10 


12.02 


-1.74 


111 


-0.03 


-0.34 


1.75 


12.01 


-1.86 


Model 


a 


P 


7 




Agal.O 


I 


-0.01 


-0.39 


2.80 


12.01 


-1.55 


11 


-0.08 


-0.42 


2.70 


12.00 


-1.48 


111 


-0.02 


-0.45 


0.95 


11.44 


-1.44 



determines how fast the transition is. The best fit values of 
these parameters are given in Table [3] 

We start our discussion with model I, which has no 
feedback (left panels). For low mass haloes the galaxy mass 
fraction m ga i is close to the universal baryon fraction, /bar — 

0. 17. This is because cooling is very efficient in low mass 
haloes. Above a halo mass of M v i r — 10 12 Mq the galaxy 
mass fraction drops significantly, due to the inefficiency of 
cooling in high mass haloes. We hereafter refer to this mass 
scale as the cooling threshold. 

The effect of feedback is to remove cold gas from the 
galaxy-halo system. The efficiency with which feedback can 
eject gas depends on both the depth of the potential well, the 
amount of star formation, and on the wind model. The net 
effect in both energy and momentum wind models is for mass 
to be lost preferentially in lower mass haloes (i.e. potential 
well dominates over star formation efficiency at fixed halo 
concentration and angular momentum parameters). Note 
the fraction of mass lost varies smoothly with virial mass, 

1. e. there is no threshold for mass loss, as for example there 
would be in a constant wind velocity model. However, the 
energy driven wind model is much more efficient at remov- 
ing baryons from haloes below M v i T ~ 1O 12 M0. This results 
in very different scalings between m ga i and m s tar with M v i r 
for the two feedback models. The parameters of these scal- 
ings are given in Table [3] In principle these differences are 
testable with galaxy-galaxy weak lensing and/or satellite 
kinematics measurements of virial masses, combined with 
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Figure 9. Dependence of mass fractions and spin parameters on virial mass for our 3 feedback models: no feedback (left); momentum 
feedback (middle); energy feedback (right). The lines show the 16th and 84th percentiles in bins of width 0.25 dex in M v ; r , the error 
bars show the Poisson error on the median. The upper panels show the mass fractions of the galaxy, m ga i = M gal /M v j r (black, solid), 
stars, m s tar = M B ta,r /M v i T (blue, long-dashed), and cold gas m co id = iWcold/-^vlr (g re en, short-dashed). The dotted lines show the 
cosmic baryon fraction, /ban as weu as /bar/2, and /bar/4- The red lines show double power-law fits to the medians, the parameters 
of these fits are given in Table [3] The lower panels show the spin parameters of the galaxy, A ga i = (i g al/ m gal)A (black, solid), stars, 
Aatar = (jstar /m s tar ) A, (blue, long-dashed) , cold gas Acoid = (jcold/ m cold)-^ (g reen > short-dashed), and dark matter, A, (black, dotted). 



measurements of stellar masses and neutral hydrogen gas 
masses. For both feedback models the maximum galaxy for- 
mation efficiency (defined as m ga i//bar) occurs around a 
virial mass of M v ir — 1O 12 M0, and is ~ 0.50 for the momen- 
tum driven wind and ~ 0.35 for the energy driven winds. 

Fig(9]also shows the cold gas mass fraction m co id, which 
is just the difference between m ga i and m st ar- Below the cool- 
ing threshold, m co id is almost independent of M v j r , while 
above the cooling threshold m co id strongly decreases with 
Mvir- These trends of m co id with M V [ T are qualitatively sim- 
ilar for all three feedback models. Thus feedback has a much 
stronger impact on the stellar mass fraction, m sta r, than on 
the cold gas mass fraction, m co id. 

4.2 Galaxy Spin Parameters 

The assumption that (j ss .i/m gsi i) = 1 underlies most ap- 
plications of the Mo, Mao, & White (1998) disk structure 
models, as well as observational attempts at measuring the 
halo spin parameter (e.g. van den Bosch, Burkert & Swaters 
2001). Thus an important question is whether this assump- 
tion is valid in galaxy formation models that include inflows 
and outflows. 

The lower panels of Fig. [9] show the spin parameters of 
our 3 models as a function of virial mass. We first discuss 
the model without feedback. For haloes below the cooling 
threshold the spin parameter of the galaxy, A ga i, is almost 



the same as that of the dark matter, which is because almost 
all of the baryons have cooled, and thus they bring in almost 
all of the angular momentum. However, for haloes above the 
cooling threshold, A ga i decreases by about a factor of 2 from 
M v ir = 10 12 to 10 13 M o . This is because the highest angular 
momentum material virializes at low redshifts, and this gas 
does not have time to cool. 

As would be expected the stellar mass fraction and spin 
parameters are lower than the corresponding parameters for 
the total galaxy. The lower stellar spin parameters are due 
to the star formation law, which causes stars to form less 
efficiently at larger radii (where the gas density is lowest). 
Thus the gas disk is more extended that the stellar disk. 

An interesting result of our models with feedback is that 
the galaxy spin parameter can be higher than the halo spin 
parameter. This is because mass can only be ejected from 
radii where there is star formation, and these are biased 
towards small radii, where the specific angular momentum 
of the gas is lower than the average. As shown below (in 
section |4~3)) , at a given halo mass, feedback is more efficient 
at removing baryons from higher surface density disks. Thus 
the increased star formation efficiency over comes the deeper 
potential well. Note that the scatter in the galaxy and stellar 
spin parameters is roughly equal to the scatter in the halo 
spin parameter, and this does not change with virial mass. 

Although the spin parameters of the galaxy and stars 
are typically different from that of the halo, Fig. [10] shows 
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Figure 10. Dependence of ratio between the specific angular momentum of the baryons/stars (black/blue) and dark matter on virial 
mass for our 3 models. The lines show the 16th and 84th percentiles in bins of width 0.25 dex in M v j r . The scatter about the median 
relation are given by <Xg a i for O'gal / m gal ) i an d ""star for (j'star/mstar)- The dotted line corresponds to the specific angular momentum of 
the baryons and stars being equal to that of the dark matter. 



that the ratio of the specific angular momentum of the 
baryons to that of the dark matter, (j ga i/m ga i), at a given 
virial mass has very little scatter (0.01 dex for model I, 0.02 
dex for model II, and 0.04 dex of model III). The small scat- 
ter holds for all halo masses, even for galaxies that have lost 
more than half of their baryons. The small scatter is related 
to the fact that there is very little scatter in the baryon 
mass fraction at a given halo mass. By contrast there is a 
much larger scatter (0.07-0.08 dex) in the ratio between the 
specific angular momentum of the stars and the dark mat- 
ter. This is due to the dependence of global star formation 
efficiency on galaxy surface density. 

Given that the assumption that (jgai/m-gai) = 1 is vi- 
olated in all of our models, care should be taken when in- 
terpreting models and observational results based on this 
assumption. However, in our models there is very little scat- 
ter in (jgai/wigai) at a fixed halo mass. This suggests that 
the scatter in the galaxy spin parameter (which is in prin- 
ciple observable), at a given halo mass, may be an accurate 
reflection of the scatter in the halo spin parameter (which is 
a prediction of ACDM, but not directly observable). 

4.3 Is Galaxy Mass Fraction Correlated with Halo 
Spin? 

Fig.|9]shows that at a given virial mass there is only a small 
scatter in galaxy or stellar mass fractions. For the momen- 
tum driven wind model the scatter in m ga i and m s tar is — 8% 
and ~ 20% respectively. For the energy driven wind model 
the scatters are ~ 20% and ~ 33% respectively. 

We now turn to the question of where this scatter comes 
from. The upper panels of Fig. Illl show the residuals of the 
m ga i — M v i r relations vs scatter in the halo spin parameter. 
For the no feedback model, the baryon fraction is deter- 
mined by the efficiency of cooling. In our model, the cooling 
efficiency is to first order determined by the halo mass. The 
very small scatter in m ga i results from scatter in the halo 
concentration, which effects the cooling in two ways. The 
halo concentration determines the MAH (i.e. low concentra- 
tion haloes collapse later, and thus there is less time for the 
baryons to cool). The halo concentration effects the density 



of the hot gas, which directly effects the cooling time. For 
the energy and momentum wind models there is a positive 
correlation (r ~ 0.6) between m ga i and A at a given M v i r . 
Thus at a given virial mass, galaxies that have lower spins, 
(i.e. higher surface density disks), are more efficient at re- 
moving baryons, despite the deeper potential well. 

The lower panels of Fig. [TTJ show the residuals of the 
wistar — M v i r relations. All models show an anti-correlation 
between the residuals, i.e. galaxies that form in lower spin 
haloes are more efficient at turning their cold baryons into 
stars. This effect is expected from the density dependent 
star formation recipe we adopt. 

These results have implications for the scatter in the 
Tully-Fisher relation. As discussed in §1, one of the surpris- 
ing aspects of the Tully-Fisher relation is that the scatter is 
independent of surface brightness, or equivalently the scat- 
ter in the VL relation is independent of the scatter in the 
RL relation. Firmani & Avila- Reese (2000), van den Bosch 
(2000), and Dutton et al. (2007) have shown that this could 
partially be explained by the dependence of gas fractions on 
surface density, such that lower surface density disks have 
higher gas fractions. Gnedin et al. (2007), on the other hand, 
invoked a correlation between disk mass fraction and disk 
surface density to reduce the correlation between the VM 
and RM relations. This works as follows. At a given stel- 
lar mass there is a range in disk sizes. Smaller disks should 
result in larger V2.2, both because the baryons contribute 
more to V2.2, and because the halo contribution increases 
due to halo contraction. However, if the smaller disks live 
in lower mass haloes, then the reduced contribution of the 
halo compensates for the increased contribution from the 
disk. This should result in a negative correlation between 
m ga i and A, opposite to what we find in our models. Gnedin 
et al. (2007) speculated that feedback would be less efficient 
in higher surface density disks, presumably because the po- 
tential well is deeper. However, in our models, the reverse is 
the case. At a given halo mass higher surface density disks 
are more efficient at removing baryons because there is more 
energy (or momentum) input from SN due to the higher star 
formation rates. 
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Figure 11. Correlation between scatter (at a given halo mass) in halo spin, with scatter (at a given halo mass) in galaxy mass fraction 
(upper panels), and stellar mass fraction (lower panels). The black lines show the median, 16th, and 84th percentiles of the models. The 
scatter in the galaxy mass fraction is correlated with halo spin (i.e. feedback is more efficient in higher density galaxies), whereas the 
scatter in stellar mass fractions is anti-correlated with halo spin (i.e. star formation is more efficient in higher density galaxies). The 
red lines shows the mean (solid) and lc scatter (dashed) of a fit of the form: Alogm = b logA/A. The slope, 6, scatter, <j (in log 10 
units), and the correlation coefficient, r, are given for each fit. The correlation coefficient is given by r = bcr x /cry, where x = log A/A and 
y = A log m\M v i r 



4.4 The Dependence of Gas Fractions on 
Feedback 

Observationally it is known that the cold gas fraction is 
higher in lower mass galaxies (e.g. McGaugh 1997 & de Blok; 
Kannappan 2004; Avila- Reese et al. 2008). Here we use data 
from Garnett 2002, who compiled B-band magnitudes, B- 
V colors, atomic gas mass and molecular gas mass for 31 
spiral galaxies and 13 irregular galaxies. We compute stellar 
masses using the following relation from Bell et al. (2003): 
log(M st ar/-M = -0.941+1.737(B-y)-0.1, where the-0.1 
corresponds to a Chabrier IMF. The gas fractions, defined 
as /gas = Af co id/(M co id + Mstar), versus stellar masses thus 
derived are plotted in green in Fig. 1121 A linear fit gives the 
following relation between the gas fraction and logarithm of 
stellar mass: 

/ gas = 0.374 - 0.162 flog ^ star - lo) (45) 

with a scatter of 0.11 in / ga s- Note that the ratio between 
atomic and molecular is a strong function of stellar mass. 
Massive galaxies have roughly equal amounts of atomic and 
molecular gas, while the molecular gas fraction is negligible 
in galaxies with M sta r ^ 10 10 Mq. Thus ignoring the molec- 
ular gas significantly underestimates the gas fractions at the 
high mass end. 



Fig. Q2] shows the gas fraction vs stellar mass relation 
for our three models at redshift z = 0. The models are given 
by the grey points. The model with no feedback (left panel) 
results in gas fractions that are too low compared to the 
observations (green points and lines). This is a result of the 
disks being too high density, which causes star formation 
to proceed too fast. Both models with feedback result in 
gas fractions that are in good agreement with observations 
(slope, zero point and scatter). In particular, the momentum 
driven wind model is almost indistinguishable from the ob- 
servations. Distinguishing between these two feedback mod- 
els would require a more complete sample of galaxies with 
robust atomic and molecular gas masses. 



5 THE DEPENDENCE OF THE MASS 

METALLICITY RELATION ON FEEDBACK 

The relation between stellar mass (or luminosity) and gas 
(or stellar) phase metallicity has long been thought to hold 
important clues to the nature of galaxy outflows. Indeed, 
the slope of the mass metallicity relation at z ~ 2 has been 
used to argue in favor of momentum driven outflows over 
energy driven outflows and no outflows (Finlator & Dave 
2008). However, the energy driven wind model adopted in 
Finlator & Dave (2008) assumed a constant mass loading 
factor, 77 = 2, and a constant wind velocity V w = 484 km s _1 . 
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Figure 12. Dependence of the gas fraction - stellar mass relation at redshift z = Oon feedback. The green points show data from Garnett 
(2002), with green lines showing the mean and 1<t scatter. The grey points show the models. The black lines show the 16th and 84th 
percentiles in bins of width 0.25 dex in M a t ar . The red lines show a fit of the form: / gas = a, + 6(log M s t ar — 10). The fits are performed 
over the range that the red lines are plotted. The zero point, slope, and scatter for each fit are given in the upper right corner of each 
panel. Note that models and data include both atomic and molecular gas. The model without feedback under predicts the gas fractions, 
whereas the models with feedback provide good matches to the observations. 
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Figure 13. Dependence of the Mass Metallicity Relation at z = 2.26 on feedback. The black lines show the 16th and 84th percentiles 
of metallicity in bins of width 0.25 dex in Mstar. The red lines show a fit of the form: \og(Z co \d/ Zq) = a + fe(logM s t ar — 10). The zero 
point, slope, and scatter for each fit are given in the lower right corner of each panel. The green points show observational data at 
z = 2.26 ± 0.17 from Erb et al. (2006). The model without feedback over predicts the metallicities, whereas the models with feedback 
provide good matches to the observations. 



As noted by Finlator & Dave (2008), this constant wind 
model is only one possible implementation of energy driven 
outflows. Thus an interesting question is whether the energy 
driven outflow model that we have implemented is able to 
reproduce the slope of the mass metallicity relation at z ~ 2. 

Fig. [T3] shows the stellar mass-metallicity relation for 
our three models at redshift z = 2.26. The metallicity 
used here is the metallicity of the cold gas, which is cal- 
culated as the sum of metals in the cold gas divided by the 
mass of cold gas. The model with no feedback has a mass- 
metallicity relation with very shallow slope of b ~ 0.06, and 
a mean metallicity close to that of the yield (i.e. Solar). The 
green points with error bars show the observational results 
for the gas phase oxygen abundance vs stellar mass from 
Erb et al. (2006), assuming a solar oxygen abundance of 
12 + log(OfH) — 8.66. Our no feedback model is inconsis- 
tent with the observations, both in terms of slope and zero 



point. The models with feedback result in steeper slopes of 
the mass-metallicity relations: b ~ 0.17 for the momentum 
wind model and b ~ 0.32 for the energy wind model. Both 
of these models are in qualitative agreement with the obser- 
vations. Overall the energy driven model provides a better 
match to the observed slope, but given the systematic uncer- 
tainties in metallicity measurements it would be premature 
to strongly distinguish between the energy and momentum 
outflow models. 



Finally we note that all of our models result in mass 
metallicity relations with small scatter ~ 0.06 dex. Thus 
we conclude that outflows are not responsible for the small 
scatter in the mass-metallicity relation. The origin of this 
small scatter, as well as the evolution of the mass metallicity 
relation will be discussed in a future paper. 
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Figure 14. Mass weighted wind velocity, Vw, versus galaxy rotation velocity, V2.21, stellar mass, M s tar, and star formation rate, M s tar, 
for energy driven (lower panels) and momentum driven (upper panels) wind models. The black and grey lines show the 16th and 84th 
percentiles of the models at 2 = 0.0 (light grey), 2 = 1.4 (dark grey) and 2 = 3.0 (black). The red lines show fits to the 2 = 1.4 models. 
The fits are of the form: log Vw = a + ft(log V2.21 — 2), logVw = a + ft(logM s tar — 10), and log Vw = a + ft(log AfsFR — !)■ The slopes 
(ft), zero points (a) and scatter (a) of these fits are given in the lower right corner of each panel. The green points with error bars show 
observational data at 2 ~ 1.4 from Weiner et al. (2009), see text for further details. Both energy and momentum driven wind models 
predict similar scalings, that are broadly in agreement with the observations. 



6 HOW DO GALAXIES LOSE THEIR GAS? 

Having established that in models with feedback, galaxies 
lose a significant fraction of their accreted baryons, we now 
investigate when and how the baryons are lost. The left pan- 
els of Fig. [T3] show the mass weighted wind velocity, V~ w , 
versus the rotation velocity at 2.2 disk scale lengths, V2.21, 
for models at redshifts, z = 0.0, 1.4, and 3.0. This relation 
has a slope ~ 1, small scatter, and holds at all redshifts. For 
both energy and momentum driven winds the mean wind 
velocity is roughly 2.5-3 times V2.2I- However, this relation 
is no surprise, since by construction we assume that the lo- 
cal wind velocity is equal to the local escape velocity. The 
non-zero scatter in the relation between Vw and V2.21 is due 
to the fact that, in our model, mass is ejected from a range 
of radii, and hence a range of escape velocities. 

The middle panels of Fig. [TJ] show the relation between 
wind velocity and stellar mass. For galaxies in the mass 
range 10 9,6 < Mstar < 10 10 ' 8 M Q this relation has a slope 
of ~ 0.33 ±0.01 for all redshifts, and both wind models. For 
the energy wind model this relation has a shallower slope at 
low stellar masses. The slope and small scatter of this rela- 
tion is a consequence of the tight correlation between VW and 
V2.21, and the small scatter in the stellar mass TF relation 
(Figs. [6]&[7]). The zero point of the relation between wind 
velocity and stellar mass evolves with redshift, which is a 
reflection of lower stellar masses at a given rotation velocity 
in higher redshift galaxies. Note that this implies evolution 
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Figure 15. Mass loading factor versus mean (outflow mass 
weighted) wind velocity for our energy and momentum driven 
wind models at redshifts 2 = 0.0, 1.4 and 3.0. The red lines are 
as in Fig. [4] 
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Figure 16. Ratio between rate of cooling (green), star formation (blue) and outflow (black) to the rate of baryon accretion as a function 
of virial mass. The upper panels are for momentum driven winds, the lower panels for energy driven winds. The panels from the left 
to right are for models at redshift, 3.0, 1.4, and 0.0, respectively. For haloes below ~ 1 X 10 12 Mq, the cooling rate is essentially equal 
to the baryon accretion rate. The relative star formation rate increases with halo mass, whereas the relative outflow rate shows exactly 
the opposite trend. Thus globally the outflow rate is dominated by the depth of the potential well, rather than the efficiency of star 
formation. Outflow rates are higher for energy driven winds due to the higher mass loading factor, this results in lower mass and surface 
density galaxies, and hence lower relative star formation rates. 



in the zero point of the stellar mass TF relation, which will 
be discussed in a future paper. 

The right panels of Fig. [14] show the relation between 
wind velocity and star formation rate. This relation has a 
slope ~ 0.38±0.01 at z — 1.4 for both wind models. Galaxies 
at z = 1.4 and z = 3.0 have similar zero points, but galaxies 
at z — 0.0 have lower star formation rates at a given wind 
velocity. 

Over-plotted in red in Fig. [14] are the observational re- 
sults at z ~ 1.4 for the wind velocity versus stellar mass 
and star formation rate from Weiner et al. (2009). These re- 
sults are based on a sample of 1406 galaxy spectra from 
the DEEP2 (DEEP Extragalactic Evolutionary Probe 2) 
redshift survey (Davis et al. 2003) which have both [O II] 
3727 emission lines (to provide secure redshifts) and Mg II 
AA 2796, 2803A absorption lines (to probe low ionization out- 
flowing gas). For the wind velocity we adopt the velocity 
where Mg II absorption is at 10 % of the continuum. 

The slope of the wind velocity stellar mass relation is 
shallower in the data than our models. However, the slope 
of the wind velocity star formation rate is consistent. The 
differences between the slopes of the wind velocity - stel- 
lar mass relation of the models and data may be caused 
by a bias against low mass low star formation rate galax- 
ies in the observations. Martin (2005) find a similar scaling 
between wind velocity and SFR in low redshift ultra lumi- 
nous infrared galaxies. Thus both energy and momentum 



wind models predict scalings that are broadly consistent 
with the observations. This agreement is primarily due to 
our assumption that the wind velocity is equal to the local 
escape velocity. It is surprising that such a naive assump- 
tion, which is most likely incorrect, is in good agreement 
with the observations. 

We have shown that the energy and momentum driven 
wind models predict very similar scalings between wind ve- 
locity and stellar mass and star formation rate. The slopes of 
these relations are only significantly different at low masses 
and low star formation rates. Fig. ll5l shows the mass loading 
factor vs wind velocity for the momentum and energy driven 
wind models. The relations from our galaxy formation mod- 
els are independent of redshift, and close to the expected re- 
lations (shown as red lines in Fig. I15|l for a model where the 
wind velocity is independent of radius. For a wind velocity 
of ~ 100 kms" 1 the mass loading factor of the energy wind 
model is an order of magnitude higher than that of the mo- 
mentum wind model. Thus observationally, it may be easier 
to distinguish between energy and momentum driven winds 
by measuring the mass loading factors of low mass galaxies. 

6.1 Relation Between Outflows, Inflows and Star 
Formation 

Fig. [16] shows the ratio between rate of (cooling, star for- 
mation, and outflow) to that of baryon accretion, versus 
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halo mass, at redshifts z = 3.0 (left), z = 1.4 (middle), and 
z — 0.0 (right). In our model the baryon accretion rate is the 
cosmic baryon fraction times the halo accretion rate, which 
is a function of the redshift zero halo mass and concentra- 
tion. Thus higher mass halos have higher rates of baryon 
accretion. 

Cooling is efficient in haloes below ~ 1O 12 M0, so that 
the rate of cooling is almost equal to the rate of baryon 
accretion. Above ~ 1O 12 M0 the cooling efficiency drops, and 
thus the rate of cooling drops significantly below the rate of 
baryon accretion. The scatter in the cooling rate is caused by 
scatter in the halo concentration, which effects the density 
of the halo and hence temperature of the hot gas. 

The star formation rate, relative to the cooling rate in- 
creases with halo mass, with a much stronger dependence in 
energy driven winds than momentum driven winds. Above 
1O 12 M0 the relative star formation rate decreases, but this is 
due to the decline in the cooling rate, rather than inefficient 
star formation. The relative outflow rate shows the opposite 
trend. Outflows are more efficient in lower mass haloes, this 
is despite the lower relative star formation rates. For haloes 
with M v i r < 10 12 Mq the energy wind model has higher out- 
flow rates than the momentum wind model. This is due to 
the higher mass loading factor in the energy wind model, 
which can be of order ~ 10 in low mass haloes. 



7 SUMMARY 

We use a disk galaxy evolution model to investigate the im- 
pact of mass outflows (a.k.a. feedback) on disk galaxy scal- 
ing relations, galaxy mass fractions and spin parameters. 
Our model follows the accretion, cooling, star formation and 
ejection of mass inside growing dark matter haloes with cos- 
mologically motivated angular momentum distributions. In 
our models the surface density profile of the baryonic disk is 
determined by the specific angular momentum distribution 
of the cooled baryons and by the assumption of centrifu- 
gal equilibrium. The surface density profile of the stellar 
disk is then determined by the efficiency of star formation 
with gas density. We model star formation with a Schmidt 
law on the dense molecular gas. We compute the molecu- 
lar fraction using the pressure based prescription in Blitz & 
Rosolowsky (2006). We consider both energy and momen- 
tum driven galaxy wind models. For both models we assume 
that the wind velocity is equal to the local escape velocity 
(from the disk and halo). This assumption maximizes the 
amount of mass loss. Our main conclusions are summarized 
as follows: 

• Veloctty-Mass-Radius Scaling Relations: Models with- 
out feedback result in disks that, at a given stellar mass, are 
too small and rotate too fast (Figs. With increasing 
feedback efficiency disks become larger, and rotate slower 
(Fig. [2j . With high feedback efficiency, disks are too large 
at a given stellar mass. However, the offset from the TF re- 
lation is almost independent of feedback, because feedback 
reduces both the stellar mass and rotation velocity (Fig. [2}. 

• Tully-Fisher Zero Point: Models with halo contraction 
over predict the zero point of the stellar mass-velocity rela- 
tion, independent of the feedback efficiency (Figs. [2]&[3]). 
Models without halo contraction result in a better agree- 



ment, but still over predict the rotation velocities at high 
stellar masses (Figs. [2] & [3| . 

• The Baryonic Tully-Fisher Relation: The Baryonic 
Tully-Fisher relation, defined as the relation between the 
(cold) baryonic mass and the circular velocity at large ra- 
dius in the gas disk is only weakly sensitive to the feedback 
model (Fig. |5J. In particular, our energy driven feedback 
model (which results in galaxy mass fractions decreasing 
with decreasing halo mass) yields an almost identical slope 
to the no feedback model (which in which galaxy mass frac- 
tions are constant with halo mass). 

• Gas fractions vs Stellar Mass: Models without feedback 
predict gas fractions that are too low (Fig. I12[) . which is a 
result of the disks having densities that are too high, which 
in turn results in star formation being too efficient. Models 
with feedback predict gas fractions in good agreement with 
observations (Fig. IT2"]| . 

• Mass Fractions vs Halo Mass: Without feedback, cool- 
ing is very efficient below a halo mass of ~ 10 12 M Q (Fig.©. 
Below this mass, the galaxy formation efficiency (defined as 
the percentage of universal baryons that end up as stars and 
cold gas, i.e. m ga i//bar) is constant at ~ 95%, while above 
this mass, galaxy formation efficiency decreases. In both en- 
ergy and momentum feedback models mass is more easily 
ejected from lower mass haloes, resulting in galaxy mass 
fractions that increase with halo mass (below the cooling 
threshold). Maximum galaxy formation efficiencies occur at 
a virial mass of M V \ T ~ 10 12 Mq. Maximum efficiencies as low 
as ~ 35% can be produced with energy driven winds with 
a feedback efficiency of 0.25. However, even with a feedback 
efficiency of 1, momentum driven winds result in maximum 
galaxy formation efficiencies of ~ 50%. See below for further 
discussion. 

• Mass Fractions vs Spin Parameter: At a given halo 
mass, higher density disks are more efficient at removing 
mass (Fig. II 1 P : i.e. energy /momentum input is more impor- 
tant than depth of potential, inconsistent with assumption of 
Gnedin et al. (2007). At a given halo mass, star formation is 
more efficient in galaxies with higher surface densities, thus 
gas fractions are lower. This helps reduce the surface density 
dependence of the stellar mass TF relation relative to the 
baryonic TF relation (as argued by Firmani & Avila-Reese 
2000; van den Bosch 2000; Dutton et al. 2007). 

• Spin Parameter of Baryons vs. Dark Matter: In haloes 
with masses lower than ~ 1O 12 M0, the galaxy spin is higher 
than halo spin, because feedback preferentially removes low 
angular momentum material (Fig. . In haloes more mas- 
sive than ~ 1O 12 M0 the galaxy spin is lower than halo spin, 
because that gas that has not cooled has the highest spe- 
cific angular momentum. At a given halo mass the spin of 
the galaxy is tightly correlated with the spin of the halo, i.e. 
the parameter (j gs ,i/m gs ,i) has a scatter of only ~ 0.02 — 0.04 
dex (Fig. I10p . This tight correlation is true even for galax- 
ies that have lost over 80% of their baryons. However, due 
to the density dependence of star formation, the parameter 
O'star/wiatar) has a larger scatter of ~ 0.08 dex. See below 
for further discussion. 

• Spin Parameter of Baryons vs. Stars: The spin of the 
stars is always less than the spin of the baryons because star 
formation is more efficient at smaller radii. Below M V [ T ~ 
1O 12 M0, the median spin parameter of the stars is roughly 
independent of halo mass (Fig. [9j . 
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• Wind Velocity vs Galaxy Observables: Both feedback 
models result in tight correlations between the mean wind 
velocity, V w , stellar mass, M s t ar , and star formation rate, 
Mst ar : Vn oc 3V 2 .2i oc M° t £ 3 oc M°iir (Fig- EE!- The first 
relation is essentially construction, since we assume the 
wind velocity is equal to the local escape velocity from the 
galaxy-halo system. The second and third relations are non- 
trivial. The scaling relations between wind velocity and stel- 
lar mass/star formation rate are broadly consistent with 
observations at z ~ 1.4 from the DEEP2 redshift survey 
(Weiner et al. 2009). 

• Differences between Energy and Momentum Driven 
Winds: The main difference is the mass loading factor, 
which scales like V^ 2 for energy driven winds, and V^f 1 
for momentum driven winds. For equal fractions of the ini- 
tial energy /momentum from the SN that drives the out- 
flow, energy driven winds have higher mass loading factors 
for all relevant wind velocities. The differences between the 
mass loading factors increase with decreasing wind velocity 
(Fig. [4j . Thus energy driven winds are much more efficient 
at removing mass from lower mass haloes (Fig. [9)l . This has 
at least three observational consequences: 1) different slopes 
of the size-stellar mass relation (~ 0.14 for energy Fig. [7] 
and ~ 0.28 for momentum, Fig. [6]); 2) different slopes of the 
relations between m ga i and m star with M V1T (see Table 2); 
and 3) different slopes in the metallicity-stellar mass rela- 
tion at z ~ 2 (~ 0.17 for momentum and ~ 0.33 for energy, 
Fig. 1 13 p . See below for further discussion. 

• Outflow vs Inflow: The median mass outflow rate rela- 
tive to the median mass inflow rate decreases with increasing 
halo mass, whereas the median star formation rate relative 
to the median inflow rate increases with increasing halo mass 
(Fig. I16p . Thus, globally, the depth of the potential is more 
important than the energy/momentum input from super- 
nova. 



7.1 Comments on Galaxy Spin vs Halo Spin 

The assumption that A ga i = A (i.e. j ga i/m ga i = 1) un- 
derlies almost all analytical and semi-analytical models of 
disk galaxy formation (e.g. Dalcanton, Spergel & Sum- 
mers 1997; Mo, Mao & White 1998; Somerville & Primack 
1999; Firmani & Avila-Reese 2000; Cole et al. 2000; Croton 
et al. 2006; Dutton et al. 2007; Somerville et al. 2008). We 
have shown that this assumption is no longer valid in models 
with outflows. A similar conclusion was reached by Mailer 
& Dekel (2002). 

The result that A ga i is significantly higher than A in low 
mass haloes helps to resolve the puzzle surrounding the spin 
parameters of bulge-less dwarf galaxies. Using observations 
of 14 late-type dwarf galaxies van den Bosch, Burkert & 
Swaters (2001) found the distribution of A ga i to have a me- 
dian of ~ 0.06. D'Onghia & Burkert (2004) measured the 
spin parameter of dark matter haloes that are most likely to 
host bulge-less disk galaxies, and found a median spin pa- 
rameter A ~ 0.028 (after correcting to the energy definition 
of halo spin). D'Onghia & Burkert (2004) assumed that in 
the best case scenario A ga i = A, and thus there is a discrep- 
ancy of a factor of ~ 2 between the observed and predicted 
spin parameters. However, with energy driven feedback, our 
models produce galaxy spin parameters a factor of 2 higher 



than the halo spin parameters in low mass haloes, thus re- 
solving the discrepancy. 

7.2 Comments on Energy vs Momentum Driven 
Winds 

In our models energy and momentum driven winds result 
in significantly different slopes of the relations between disk 
size and stellar mass. Observations of the size-stellar mass 
relation for late-type galaxies from Shen et al. (2003) find 
a slope of 0.14 at the low mass end, which favors our en- 
ergy wind model over our momentum wind model. How- 
ever, at the high mass end, a number of authors find steeper 
slopes (e.g. Shen et al. 2003; Pizagno et al. 2005; Courteau 
et al. 2007), which is in better agreement with our momen- 
tum wind model. However, there are a number of uncertain- 
ties in the observations (such as determinations of stellar 
masses, inclination effects on galaxy sizes) which need to be 
quantified before firmer conclusions can be made. 

Our energy and momentum driven wind models also re- 
sult in significantly different slopes to the relations between 
galaxy mass, stellar mass and halo mass. These relations 
can in principle be directly tested with galaxy-galaxy weak 
lensing and/or satellite kinematics. 

7.3 Comments on Why Galaxy Formation is 
Inefficient 

Observations of halo masses from weak lensing studies (e.g. 
Hoekstra et al. 2005; Mandelbaum et al. 2006) and meth- 
ods that match the stellar mass function to the halo mass 
function (e.g. Yang et al. 2007; Conroy & Wechsler 2008) 
find that the peak galaxy formation efficiency has to be rel- 
atively low ~ 0.33. We have shown that low galaxy forma- 
tion efficiencies are also required to explain the zero points 
of the relations between velocity, stellar mass, and size of 
disk galaxies. 

We have shown that mass ejection through supernova 
driven winds provides at least a partial explanation for this. 
However, by assuming that the wind velocity is equal to the 
local escape velocity, results in maximal outflow rates for 
a given energy /momentum input. While such a scaling of 
wind velocity with galaxy escape velocity is at least sup- 
ported by observations at low and high redshift, galactic 
winds will likely have a range of velocities, which will re- 
duce the outflow rates from those in our model. Needless to 
say, the scaling between wind velocity and escape velocity, 
as well as the mass loading factor, need to be investigated 
further with hydrodynamical simulations. 

Even though we have adopted a maximally efficient 
mass outflow model, we still need to use 25 percent of all 
SN energy (or 100 percent of all SN momentum) in order 
to eject enough mass to fit the data. It remains to be seen 
whether such high efficiencies are realistic or not. This re- 
quires detailed hydro-dynamical simulations with radiative 
transfer, that accurately model the complicated multi-phase 
structure of the ISM. It is likely that one needs to invoke ad- 
ditional mechanisms to explain the low baryonic mass frac- 
tions observed in galaxy mass haloes. 

An alternative explanation for low galaxy formation effi- 
ciencies is that most of the baryons never accrete onto galax- 
ies in the first place. In massive haloes M v i r > 1O 12 M0 (in 
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which gas is heated by an accretion shock), accretion onto 
the galaxy could be suppressed with multi-phase cooling 
(Mailer & Bullock 2004) or additional heating such as from 
AGN, or mergers. However, in low mass haloes, most of the 
baryons are accreted in cold streams. Rather than disrupting 
these streams from the outside, such as with feedback from 
the central galaxy, a more likely scenario would be to disrupt 
them from the inside, i.e. by re-heating baryons and eject- 
ing them from the cold flow into the inter galactic medium, 
before the cold flow reaches the halo. We have shown that 
even though the star formation efficiency is much lower in 
lower mass haloes, the mass loading factor is typically high, 
especially for energy driven winds. Thus even small amounts 
of star formation in low mass haloes could be sufficient to 
significantly reduce the baryon accretion rate, and hence 
baryon mass fraction of the main galaxy. Coupled with out- 
flows from the main galaxy, this could result in galaxy mass 
fractions in better agreement with observations, but with 
a lower (and more realistic) conversion efficiency of SN en- 
ergy/momentum into galactic winds than required by our 
current models. 
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